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1. INTRODUCTION 


1,1 j f ne ^ al - Wit h Che maturation of the space shuttle concent 
modeq bl f launch . vehicle fo ^ earth orbital missions, two diver- 
" * 0 P erat ^ on ha ^e been defined. One mode involves the use 

packages in orbit and 08 * 3 ^ 5 VG '‘ iCle plaCing free flyin S experiment 
* c 5 ® ° b d re P lacin g, repairing or servicing existing pack- 

ages In addition it will perform a crucial role as a manned expS- 
base, remaining in earth orbit from 7 to 30 days performing 
various experiments with equipment mounted in the payload bay 

in lo^earth orbitT^ ° f ^ expGrimen ts Proposed for operation 
ir low earth orbit m conjunction with a manned vehicle indicate 

that nearly 45 percent of the payloads require pointing accuracy 

greater than that afforded by the orbiter capability using the Reac- 

the capaM? r ySte f ^ ^ 39 percent have requirements exceeding 

Moment P S ^ ^ ^ the dl ™ ct Ga ^rol of Control 

hi'"w G/roscopes (CMGs). One concept under study to meet these 
higher accuracy pointing requirements is to mount the experiment 

bi?er S cb f Pallet s *-ructure which in turn is attached to the or- 
d- r .b thrDUgh a . suspension system which isolates the pallet from the 

st^nr^pL^Kn? 8 * ^ ° rblter WhlCh P «^* oft- 


The study described in this report involves the definition and 

attached°t 0 °the ' ^ C °** to1 system ^ilizing CMC actuators 

attached to the pallet and controlling it directly, affecting and 

affected by the orbiter only through the suspension system. 

• b ’ 2 . Sl:a dy Objectives - The objectives of the suspended pallet 

pointing performance study were the following: 


a. Selection of a CMG actuator and definition of a 
figuration that best satisfies the overall mission 


system con- 
requirements . 


b. Development of a detailed nonlinear model of 
CMG for analog computer analysis of the actuator 


the selected 
characteristics . 


turban™ V£hlCl£ eontro1 law translating vehicle dis- 

ce effects into a correcting torque command and a CMG con- 

troi law to generate appropriate gimbal rate commands to generate 
the commanded torque. b ieL,iL - e 


d. Development of a model of the rigid body dynamics of 
let and orbiter connected through the suspension system. 


the pal- 
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usad I ln e r^' :i ° n i 0f aU " : ° dGls into a hybrid simulation to be 

by Cht “ P ° lnCla Z Performance realisable 

y cue cmo pallet-suspension-orbiter system. 

All of the above objectives were met during the course of the study 
^ a if 6 behavior of CMC systems and the software required 

ri.a f r ° l c s rir? r an f e su ; h a BysceK - <*?<*- 

formed 23 T** T* tha “Ven years that tSe s^dlS 1“ ?“ 

formed, and was directly applicable to the performance of the present 

Data generated during the course of the Research an H j n i • 
tions Module (RAM) study (May 1971- June 1972) in which todix ias" 

in the^area^f 11 ’ 1 '' 13 ' P ^° VldlacI villuabl ° “put to tile present study 
the area of proposed experiment payloads for the shuttle vehicle 

?n add h3 ti™ U th ementS Pl * Ced UP “ bbal vehicle La J “C P lyllt£ 
candidate^ * „!T“ Sal “ d *» tha “aalysis and desig* of 
applicable to the preset “ud^ ““ applidati °"= is Erectly 

1971-Sehl973rL^“h ^"wif? priff^tL^^ 

Vided valuable inputs to the present study LX LTol tllJT 
in 'thi e desin lr T ntS d ,“? “ ?r0i ’ OSecl astronomy experimentation and 
reaction control X lbl % a “ it "'?* s y“ta»s which included 

be used to meet rtf ’ • Hlperl '"'“ lt Isolation systems that could 

ue usea to meet these requirements . 

c^v stud y to Define Logic Associated With CMGs to Maneuver and 
Stabilize an Orbiting Vehicle" was performed by Bendix 11^9 
September 1973) defining a CMC system capable of con Collin e thl 
orbiter and developing the software required to manage and control 
the system. The experience gained in the performed of ImTII ] 
was directly applicable to the present studyf p^Sculafly in th? 
areas of CMC system selection and CMG control software. 
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nr norim! T* 1 ’*™.*™* ined during the performance of the preceedinc 

=EH“sH" : s £=saa 

asSKSSS^ff 

“ “lectio of tha CMC actuator was the avoUabiStv^rf”, , 
hardware in order to ohi-aHn cr . availability of developed 

~ = tt~&£=£s^sr 

law but also baciuL lt anlhi°fl tf SeleC 5 io " of a di Sital CMC control 

til. iu^fcontror'Jsta^aS"^ ^"^'^^‘"^^^“^actuat^s and 
quantization in the ^rf^r^al™'ti„“? ti0 " ^ ““ ° £ 

Of thfLtle^pLroflq-'f “ mi , d - 197 ' i al »»E with pallet designs 

n' J h , e ., CMG f ySCein is sized Co al l°w continuous pointing for 
one orbit under the worsi' pnnn ythd ^ • . 8 

accumulations. gravi y gradient momentum 

p' th * C , MG s y stem selected consisted of four double gimbal 

unlimited ^gimbal TrlLZ^ ^ ^ *> a11 ™ 

£ «r^ 5“ r 

an optimal, pseudo-inverse law giving decoupled vehicle control 
while minimizing CMG gimbal rates. venicie control 

c^pUa^srlriSon^nil *“*■«•« ^nal 

*--• , , uon .revels, dead zones and servo nonlinfnH- 

actuator. &V& t0 evaluate actual characteristics of the 
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1 ' 6 Limitations - This study was geared to the determina- 

tion ot a bound on the pointing accuracy of the floated pallet con- 
cept using a CMG control system. Although the pointing stability 
achievable by the CMG system selected has been determined, the point- 
ing performance that could ultimately be obtained would require a 
more comprehensive simulation for its determination. This simulation 
c-ould include vehicle flexibility, CMG vibration characteristics, and 
adequate representation of disturbances of all types. A more accurate 
determination of the pointing stability achievable by the CMG system 
is desirable since it would have a very strong effect on the type 
and complexity of the required experiment integration equipment and 
pointing devices. This in turn can appreciably effect the overall 
cost of the shuttle experiment program. 


Although tra< ? e studies were performed to select a candidate CMG 
system configuration from several alternatives, they did not consti- 
tute a major portion of the study and hence were cop level in 
nature. When the design of an actual CMG system for the pallet is 
contemplated, considerably more effort would be required in trading 
various candidate CMG system configurations in order to determine the 
optimum configuration that would meet overall pallet/orbiter require- 
ments in a cost effective manner. Mass characteristics should be up- 
dated to reflect the present configurations. Additionally the latest 
experiment payloads slated to fly in a sortie mode should be examined 
in order to better ascertain and update the requirements placed upon 
tne vehicle by these experiments . 


Finally, software required to implement the vehicle control sys 
tern was not considered beyond the point of specifying basic control 
laws. Software required to accomplish momentum management in par- 
ticular should be considered in the selection ox a CMG actuator 
cluster based on detailed mission requirements, since the momentum 
storage capability essentially determines the number of actuators 
By using knowledge of vehicle attitudes together with biasing of the 
actuator momentum states, the possibility exists of accomplishing 
some missions with fewer actuators than specified here, allowing a 
corresponding increase in payload weight. 

1,7 Suggested Additional Effort - The following additional 
investigations affecting the floated pallet control system are 
suggested : 


a. Review the instruments and payload experiment combinations 
which are candidates for the floated pallet to revise and up- 
date mission requirements as they affect the total control 
loop . 
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f lexihT^K "i ° n °5 th<2 pointiT1 8 accuracy realizable when 

J“ ble b ° dy dynam±GS > CMC vibration propagation and non- 
deal tensors are included in the analysis. This would affect 
Lhe design and specification of instrument pointing devices 
independent of the primary control loop. 

sensor^th Gandld ^ te CMG control system configurations including 

indepth P ‘ Ust ” is5io1 ' tequiremeats. Perform 

indepth trade studies between the various CMG control system 

tiofthat^iTl “ ! rder ^ determine the optimum CMG configura- 
pff “ b Wl11 meet GVera11 mission requirements in a cost 
effective manner. These studies should include detailed hard- 

SaSit'1 f d W±t J reSP£ f " "***• P-er cons^ptL, 

i G °f as wel1 as the '^ct various QIC configu- 

In additi 1 ^ 6 ° as ° ftware G ^plexity and overall software cost. 

rate control C la^ ^ .° f COntro1 lo 8 ia (e.g., CMG gimbal 

rate control laws, singularity avoidance laws, etc.) required 

to satisfactory manage the CMG control system should be eval- 

comnf t^ rGm the . Staadp ° int Gf softw are complexity and required 
computer capacity for implementation as well as overall^er- 

formance. These software trades will be of key significance 
in determining whether a dedicated computer or the centrS or- 

coSgSir W ° Uld ^ US£d ^ dr±Ve CMC system 
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2 . 


SUMMARY 


The floated pallet definition study consisted of three related 
bat separate studies. The first was an evaluation of alternate 

suits X of which m a ° rbiter experin,enl; al instrumentation, the re- 
ults of which are presented in volume I of this final report The 

second was the pallet pointing performance study treated in this 
volume II The third was a hardware conceptua/design study in! 
PM p dlnS pallet; suspension and retention systems, experimenterection 

reporr^Surfs'l’ ^ J eSultS presetlted in volume III of this 

durinv* th! 15 a r dla § rara of the logical flow of activity 

during the performance of the pointing performance study. 

The basic requirements under which this study was performed are 
own as block 0 on the diagram. The contract statement of work 

V aSkS and I-put data for the studfi”! 

eluded baseline shuttle data, particularly payload interface in- 

daSTinc! d raqUlred f ° r Pallet and mass properties. Pallet 

tvnic!! i - maSS propertles » structural and dimensional data and 
this =(-. xp ? ri ™ eat layouts. The mission requirements applicable to 
this study included orbital parameters, vehicle attitudes and typical 

SStXnf COnStrallLtS ; The Panting requirements included stated 
pointing accuracies, stability rates and associated parameters, con- 

files ° 01 i Gia ^f cp£risties desired and some typical disturbance pro- 
files under which control loop operation was to be Investigated 
This material is discussed in sections 3 and 4 of thilvolSme 

rur poxntXn § performance study was initiated with block 1 the 

Telit?, Trad£ StUdieS - The basic requirements a^d ve- 

hicle mass properties were translated into CMC system requirements 

such as torque levels, momentum storage, redundancy and actuator 
bandwidths needed to meet the overall study requirements wiS results 
g ven m section 5. Next, the industry was surveyed for available or 
nearly available CMG actuators which could meet the cluster reemire 
mencs. The candidate CMGs ware traded with respect to s“e weiSt' 
power, and cost on an actuator and cluster basis with both single ’ 
and double gimbal actuators considered. From this survey, the^iG 
actuator and cluster configuration to be used throughout the stSiy 

r^s f "w ms portion ° f *• is * -cti Ti y 

In parallel with the hardware studies presented in volume III 
control system study was performed and shown as blocks 2 3 and 

4. Initially a complex nonlinear CMG actuator model including 11- 

compliances, frictional characteristics, deadbands, and shock- 
mount characteristics was derived and programmed for analog Simula- 
tiop i,: order Co evaluate the acceptor charao ceriaties , I^paSSlar 
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Figure. 2 1. Logical Flow of the Pointing Performance Evalnation 






“r airsSifit:„r del r s si " piified ~ 

material is given in section 7 ! iL llTtltl ll ™ ™ 

initiated with ' th^def ini?LrS t±n ^ Perf ° rmanCe StUdy ° f bl ° ck 3 was 

=L ri8 h id 

derived eLL-eraor^ 1 ^’/ 116 nonlinear shoefcmoumted CMG models 

H^^SSiSSrEs 

the e- S 3 unctlon Gf CMG parameter variation. The results of 
of thi s n volume? r GrmanCe Study are shown and discussed in section 10 
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3. 


VEHICLE DESCRIPTION 


S e S’«s sr al c “ k ^ p-SSLt be 

(SRBs) f t r“,“ E1I ' S ' and tu ° Solid Pocket Boosters 

Thf&s and 5““ ?““ le Fll8ht *■> ohow, in figure 3 ”. 

thru2ffor?if^Ie bl “L^» enEl ” 6S “» *" Proving 

^ jss^aS? ‘tt'tsisfsszz ” here the “ ternai 

(OMS, is fired to piece the orbiL in the desi«d“IrMt? SUbSyStem 

- s : s LI- ^cS^: slE “ -» - -- 

1» th. 3 1e y ff^ i fMiy^ ' K; T “fj- oy the orbiter 

lsss/is ^et°»r:Sit° che ° rbiter 

be rigid bodies. Figure' 3-2 is f* considered to 

combined configuration “ U Vle ” °* «“ orbiter/pellet 

isolator suspension point is assumed to exhibit linear 

- “ a -SiS isolatioiLwith 

in tabS S 3!l 0Pe !L‘' the orbiter /pallet configuration are given 

tL InertiaL (IWoord ; '” aSS loCMi °” s »« given referenced to 
at r. ^ ^ rdrnate system, which has its origin located at 

£: orbiter ^ithles^t to^he S*"V“ 3 8hoW the of 

°h the h Payl ° ad bay iS ^“°^cSes te above^he inertiS’axis 

SSLlli, ^ ° f 811 T' 

attitnd!' 1 f rb ^ ter " ® nce the orbiter has reached the desired orbit 
attrtude control as maintained by the Reaction Control System fR^ ’ 

or by some control mechanism which has been olared h ^ ^ , 

bay. Through usa of the RCS system, 
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Figure 3-2. Orbiter /Pallet Configuration 



Table 3-1. Orbiter /Pallet Vehicle Mass Properties 


Mass = 82,394 kg 

Center— of-Mass wrt I 
(27. 97 , 0. 0027, 9. 54) T 

Coordinate System ' 
meters 

I xx =l,010,359 kg-m 2 

I xy =~9,490.8 kg-m 2 

V" 7 ’ 400,759 kg_ra2 

I xz =266 » 419 kg-m 2 

^ zz = 7,614, 979 kg-m 2 

Iy 2 = “ 3 » H8 kg-m 2 


Al x“ I zz _I yy _214 ’ 220 k §-™ 2 

Al y" I xx“ I zz =_6>604 ’ 620 k g" m2 
Al z =I yy ~ I xx = 6,390,400 kg-m 2 
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LOCATION OF PALLET CM wrt ORBITER CM 

■ — v- — ' 



LOCATION OF ORBITER CM wrt I SYSTEM 


Figure 3-4. Coordinate System Location 



achieving and maintaining any attitude to a certain degree of accu- 
racy. If more accurate pointing is required for a particular pay- 
Porting capability must be included within the payload, 

. being the case for this particular study. The method chosen 
to increase pointing capability for the present study is to control 
the pointing of the pallet via use of the CMG cluster. 

rc\ P r °P erties for the orbiter are defined about the orbiter 

( ) coordinate system and are given in table 3-2, Figure 3-4 shows 
the orbiter center-of-mass location with respect to the I system. 

n 3 ‘ 2 ' 2 - P - allet - ” Experiment equipment will be mounted on the 
pal_et. When performing experiments, it is undesirable to allow 
disturbances caused by orbiter RCS firings. Therefore, the pallet 
wiU be controHed by the CMG cluster which provides pallet pointing 
*° W ^' lln ° ne S6C * S±nce the orbiter will also be controlled by 
the CMG s via the isolators, the RCS deadband should not be exceeded 
and no firings should occur. 

The pallet will consist of mechanically coupled segments and 
the paile t attached to the orbiter through the isolators. There 
will be a separate retention system which will support the pallet 
unng launch. Once in orbit, the retention system will be dis- 
engaged and the pallet will float on the isolators. 

Figure 3-5 is a pictorial diagram of the pallet. A detailed 
description of the pallet along with a detailed analysis of suspen- 
sion systems is given in volume III of this final report. 

Mass properties for the pallet are defined in table 3-3 with 
respect to the P system and the location of the pallet center -of - 
mass is shown in figure 3-4. 

3,3 QM? Control System - The CMGs are required to provide the 
pointing contro! of the pallet and to absorb the momentum buildup 
cue to GG torque. The CMGs are mounted on the pallet through iso- 
lators which provide a bandpass of 20 Hz with damping of 0.15. The 
purpose of these isolators is to prevent high frequency disturbances 
generated by the CMGs from reaching the pallet. Volume III gives a 
detailed discussion of the CMG mounting system, while the CMG con- 
trol system is discussed in other sections of this document. 

3-4 .Coo rdinate System Definition - All vehicle coordinate sv— 
tems are defined nominally with the X axis directed to the rear of 
the orbiter, the Z axis pointing out of the payload bay and the Y 
axis completing the right-handed system. 
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Table 3-2. Orbiter Mass Properties 


Mass = 71,419 kg 


Center-of-Mass wrt I Coordinate System 
| ( 28. 24, 0. 005, 9. 46) 1 meters 


J »s' 985 - 683 k E-" 2 

J yys =7, 219,756 kg_m2 

j 2zs =7,3S6 > 523 kg-m 2 


J xys = ~ 8 ’ 135 kg-m 2 
J XZS =256 ’ 250 k g-tn 2 

J yzs = ~ 4 ’ 067 k g-m 2 


Table 3-3 . Pallet Mass Properties 
Mass = 10,975 kg 

Center-of-Mass wrt I Coordinate System 
(26.22,0.01,10.05) meters 

J xxp^ 21 ’ 293 k S~ m2 J X yp=-99i kg-m 2 

J yyp =138 ’ 821 kg-m2 J xzp =-l,231 kg-m 2 

J zzp =133 > 426 k g-m 2 J yZ p =903 kg-m 2 
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Figure 3-5. Pallet System 
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3.4.1 Inertial Pointing System (X T ,Y„,Z T ) - This system is 

based on the physical system by which all shuttle components are 
defined. The origin of this coordinate system is fixed at the nose 
of the disposable boost tank. All coordinate systems used in this 
study are referenced to this system either directly or indirectly. 
It is used in this study to define the inertially held attitude. 

The Xj axis is 10.16 m (400 in.) below the shuttle payload bay 

centerline (figure 3-3) and the I system origin is 6.045 m (238 
m.) ahead of the orbiter nose. 


3.4.2 Orbiter System ( X s» Y s » Z g) “ This system origin is lo- 
cated at the orbiter center-of-mass. The origin of the system is 
translated from the origin of the inertial system by the vector 

p S * E <5 uati on 3-1 gives the rotational transformation from the I 

coordinate system to the S coordinate system through the angle t . 

S 


r t i= 


i 

-0 


Sz 


0 


Sy 


0 

1 

-0 


Sz 


-0 


Sy 


Sx 


Sx 


(3-1) 


3.4.3 Pallet System ( X p» X p> Z p) - The origin of the pallet 

coordinate system is located at the center-of-mass of the pallet. 
The pallet system is translated from the I coordinate system by 

the vector p p . The transformation from the I to the P coordinate 

system is given by equation 3-2 through the angle 0 . 




1 

- 0 . 


Pz 


0 


Py 


0 . 

1 

- 9 . 


Pz 


PX 


-0 


Py 


Px 


(3-2) 


3.4,4 P —System — Since the CMGs will not be aligned with one 
another, and the torque output from the CMG/shockmount is defined 
along the CMG base coordinate system, a set of transformations is 
required to relate the individual torque outputs to the pallet (P) 
coordinate system. A P’ coordinate system is defined to make the 
necessary transformations. Each of the P* systems is centered 
at the center-of-mass of one of the CMGs. The relationship be- 
tween the pallet (P) and P* system is discussed further in section 6 


. 3 *^* 5 G S y s tem * Identical CHGs are used in the cluster, 
.he derivation of the CMG model is given in detail in section 7. 
All coordinate systems associated with the CMG are defined in sec- 
tion 7 and the torque output from the CMG/shockmount system is de- 
fined in the P' coordinate system. 


t 
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4. MISSION REQUIREMENTS 

In this section typical floated pallet mission requirements 
are defxned that influence the conceptual design of the CMC atti- - • 
tude control system. These requirements are utilized in succeeding \ 

sections of this report to size and configure the CMC actuator ^ 

cluster and to accomplish the preliminary design of the control 
system. The data discussed in this section is based in part on 
requirements of the dedicated solar sortie and the combined astro- 
nomy missions projected for the floated pallet concept. The 
alance of the information is baseline data extracted from the 

contract statement of work and sources related to the statement 
of work. “ 

_ Vehi cle Orbit - For the purposes of this study a cir- 

cular orbit of 5Q0 km (270 n.mi.) will be used. The orientation 
ot the orbit in space does not impact this study in any way. 

.Vehicle Atti tudes - Both sortie mission payloads con- 
sidered for the floated pallet require inertial or nearly inertial 
attitudes which allow individual instruments to be directed toward 
stellar targets or toward the sun. The vehicle attitudes considered 
for this study are: 1) an inertial attitude with the composite 

vehicle X axis perpendicular to the orbital plane (XPOP} and 2) 
an inertial attitude with the X axis constrained to lie in the 
orbital plane (XIQP) . 

The astronomy and solar physics payloads can be pointed and 
stabilized using either of the above attitudes. These two inertial 
attitudes impact with the designs of the CHG attitude control and 
experiment pointing systems in different ways. For an inertial 
XPOP attitude, the vehicle's Y and Z axes are constrained to the 
orbital plane thus reducing the spacecraft's rotational degrees-of- 
freedom from three to one. The vehicle can be rotated about only 
the X axis. In order to permit a hemispherical experiment coverage 
capability at least one wide angle glmbal must be added to the ex- 
perimental instrument to provide the second degree-of -freedom re- 
quired. 

For an inertial XI0P attitude, the vehicle's X axis is con- 
s ^ ra ^ ne< ^ to t * ie orbital plane thus reducing the vehicle's rota- 
tional degrees-of-freedom from three to two. The vehicle can be 
rotated about the X axis and an axis normal to the orbital plane. 

Because the vehicle still has two rotational degrees-of-freedom, ' 
the experiment package can be pointed anywhere in the celestial 
sphere hardmounted to the spacecraft, i.e., no wide angle gimbal 
is required. The cost of this additional rotational degree-of- 
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freedom for an XIOP attitude is a 

Cum storage requirement per orbit over ^ rease . ln t h e momen- 

attitude Thi <= p c over Chat required by an XPOP 

impacts the size of 7 St ° raee requirement significantly 

related to additional CMGs requSd anfmorel * direGtly 
desaturation "dumps." Depending on the fret iuent CMG momentum 

Pointing and stabilization ^reJuiremenL Payl ° ad and its 

attitudes may be required or Llvih^ ’ ? U “ Ehesa twG iner tial 

are particularly attractLe G ?? aln characteristics that 

sions. attractive for a specific mission or class of mis- 

4*3 Vehicle Man euver Guideline d - ,-u 

f haS beeT^ififid tkat~the~CMG~c^ntrol ^ 

for any gross vehicle maneuvers Larger- a oi t 111 E be rec I u;Lre d 
accomplished (if required) bv the n gle maneuvers will be 

polling h re3 r Ct “ that «U*. «» Pallet 

of the Z axis with respect to consi erin S the rotational displacement 

pointing m „'£ ±n * I6ial "tectinn. The 

p«iietion ffgioo P °T^^ xrrrdb w de ff n r 1 f ° r the 

pallet loop consists of -t Z v 3 dl * bandw ldth of 2 Hz. The 

laws, the CMG actuator 
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5. 


CMG SYSTEM REQUIREMENTS 


. clig t orque and momentum storage requirements are defined f 
two shuttle pointing orientations, i.e., vehicle X axis perpen- 
dicular to the orbital plane (XPOP) and vehicle X axis in the 

nriate 1 / 13 ^ (XIG ^ * These ^uirements are based on the appro- 
gradlent disturbances with consideration given to 
the effects of aerodynamics, internal motion, and venting dis- 
turbance torques. Consideration is also given to the effects of 
vehicle positioning misalignments, principal axis misalignments 
and gravity gradient desaturation maneuvering on the torque and’ 
momentum requirements.. 

5,1 Obiter /Pallet Mass Properties - The orbiter and pallet 

M DE f 6d “ th±S StUdy are surnmar ized in tables 

l 3 ° f hlS V ° lume * The Properties are defined with 

respect to axis systems shown in figure 3-4. Using the basic mass 

h? f th£ tW ° GOm P° neili:s > tbe center of mass of the shuttle/ 
call ^ a ° mblnat:L ° n can be determined, then the inertia tensor of the 
combination can be computed. The mass properties of the composite 
vehicle are given in table 3-1. 1 

„ SlaCa the products of inertia are small compared to the mo- 
ments_ of inertia about the vehicle axis, the misalignment of the 
principa! axes can be estimated quite closely by assuming that 
the principal axis inertia tensor has the same values as the main 
diagonal of the vehicle* axis Inertia tensor and further aas™“g 

that a small angle rotation about the vehicle 

axes will define the transformation from vehicle to principal axis 
coordinates. s 


t T ] = 
P v 


1 

-e 


zp 


-e 


zp 


1 

-e 


yp 


“xp 


yp xp 

Then, the inertia tensor in vehicle coordinates [J ] can be 


written 


as : 


[ J V J = [ Tldiafid ,1 ,1 }i[ t 1 

V V p Try* ~ 


XX' yy ' zz L p v J 
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Performing the similarity transformation of the above and soLvine 

i-U „ „ ~b 


for the components of e" then gives: 


£ xp“ A l~ =O - 0145 57 rad=0.834‘ 
x 


(5-1) 


e vp Jj~ =0.040338 rad=2 . 31° 


(5-2) 


£ zp“ A? =0 * 0 O1485 rad=0 . 0851° 


(5-3) 


The composite vehicle mass properties to be used in the CMC sizinv 
can be taken directly from table 3-1 and equations 5-1, 5-2 and 

are * defined ^ 1 * ^ the ^ ^ fqultions 

are defined as rotations about the vehicle axes. 

aetiJ’on M T entUm Re ^ ulreTnenf ~' q " The disturbances 

?n ^ ? h ? Vehlcle whlch can be considered in the CMC sizing 
include gravity gradient, aerodynamic, internal motion, and ve- 
hicle venting torques, for the purposes of this report vehicle 

Jb? ^ g Wi * X n f C be COnsidered due t0 lack of definition. In 
addition, it should be noted that analyses of the Skylab missions 

Sfect the CMrT 1Cal VentinS dl ®turbances do not significantly 
affect the CMG torque requirements and are rarely a problem in 

momentum storage (reference 1) . Prior studies of shuttle systems 
(reference 2) indicate that peak aerodynamic torques are an order 

, Smal , ler / han the peak gravity gradient toques for 
P attitudes and of a similar magnitude for XPOP attitudes 

factoS^nlLff 8 ^ 11 thUS ^ a PP r -^ed by the ap^iate 
sWr a ? p ^ led to the gravity gradient torques for each attitude 
Period internal mass motion disturbances will not affect the 
cluster sizing, although they are considered In detail as they affect 
vehicle pointing in subsequent analyses. 7 eGt 

The gravity gradient torque acting on a vehicle is given by: 

r _ 


a a AI 
y z x 

9 

P = 3w a a AI 
gg o x z y 

a a AI 

x y z. 


(5-4) 
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Where the a^, i-x,y,z are the components of the unit vector along 
the local vertical resolved into vehicle coordinates, the AI are 

defined in table 3-1, and the products of inertia are ignored. 

This is equivalent to stating the sizing will be performed in the 
principal axis system. The basic sizing will be accomplished with 
this assumption, then the effects of angular offsets with respect 
to the basic attitudes will be developed to investigate the effect 
of principal axis misalignments and pointing offsets. In all cases 
the circular orbit of 500 km will be used, for which the orbital 
rate is: 


U G “1- 107x10 ^ rad /sec (5-5) 

Another significant contribution to torque and momentum storage 
requirements is vehicle maneuvering. Although any gross vehicle 
maneuvers will be accomplished without the use of the CMG system, 
these maneuvers will be examined as to their effect on torque and 
momentum requirements. A minimum desaturation period of one-third 
orbit will be considered with a small angle desaturation scheme 
assumed . 

5 -2.1 X Axis Perpendicular to the Orbital Plane (XPOP) - The 
XPOP orientation is advantageous in the sense of minimizing torque 
levels and accumulated momentum, however, to achieve full pointing 
capability the experimental payload must have a wide angle degree 
of freedom about one of the axes constrained to the orbital plane. 

Referring to figure 5-1, it can be seen that for XPOP, the 
components of the unit vector along the local vertical are: 


- 


* • 

a 


o 

X 



a 


COSW t 

y 


0 

a 


sinia t 

* z _ 

. 

L o J 
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The corresponding gravity gradient torque is obtained from equations 


5—4 and 5—6 as : 


t - 3 2 

I - — OJ 

eg 2 o 


AI sin2to t 
x o 

0 

0 


(5-7) 
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torque 


Substituting the 


occurs for sin 2 V - ± 1, i.n, n=?,3,5?7. _ 

numerical values for « q and 41^ gives the values of peak torques as 

+P-4 n-m at these points for the ideal XPOP atHf-iide ^ 

accumulation as a function of time is the time Integral of equation 


-*• 3 

H - f a) 
gg 4 o 


AI (l-cos2u> t) 

X o 

o 

o 
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f“roftSlr 5 "?. Sh0W , S that Che momentum accumulation as 

desaturation^would bl required" thuS 

occur** w * required . The peak momentum accumulation 

o C “ ’ “o’ 1 or 3ir/2 at which point the magnitude 

of the accumulation is 355 n-m-sec. 

The effect of offsets of the vehicle from the ideal attitude 
can be determined by defining a small rotation a. about each ve- 

verticaf^eome? C0mp0neritS of the unit vector along the local 


r = 
o 


— £ 


— £ 


e cost) t— £ sinw t 
z o y o 

costa t+£ sinoj t 
ox o 

sintu t— £ costa t 
ox o 
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The gravity gradient torque including the offset rotations is thus: 


£ 0 -v 2 

T =T -3(0 £ 
Hg gg o x 


AI cos2co t 
x o 

0 

0 


3 2 

- ~ t) £ 

2 o y 


, 3 2 

+ -X t) £ 

2 o z 


(l-cos2tj t)Al 
o y 

Al sin2co t 

Z o 


Al sin2(o t 
y o 

(l+cos2t) t)Al 
o z 


(5-10) 
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For rotations £ x there is no effect on the torque magnitude, merely 
a shift of the values in the orbit by the angle as can be seen 


from figure 5-1 and verified from equation 5-10 with e =e =0. For 

Y offsets (£^=£^=0) the peak torques increase from the^deal values 

and shift in pairs from the quarter orbit points toward u^v/2 and 

3ir/2 Figure 5-2 shows the torque magnitude history for one orbit 
for the ideal attitude and for e y -l°, 2°, and 3°. For e =1° the 

torque peaks occur at 53.6°, 126.4°, 233.6°, and 306.4° with - 

^ 8 M° U fnd°270° 5 ^ 2 ,“ ^ 3 ° the t0rqUe Peaks have fleeced 

Vnr 7 l 7 u h “^mtudes of G *85 n-m and 1.27 n-m, respectively. 

or Z offsets the results are similar with the peaks moving toward 
the nir points of the orbit. For ^=1° the peaks occur at 36.1°, 

X4 j‘o°’ and 323.9° with a magnitude of .5 n-m and for 2° 

and 3 the peaks occur at 0°, 180°, and 360° with magnitudes of 
U .03 n-m and 1.24 n-m, respectively. 

The momentum accumulation can be found by taking the time 
integral of equation 5-10: 


0 3 

H =H - - 10 e 
gg gg 2 o x 


Al sin2io t 
x o 

0 

0 


3 , 

- — 03 e 

4 o y 


L 3 

+ 7 - o) e 
4 o z 


AI (2oj t-sin 2 o 3 t) 
y o 0 ' 

Al (l-cos2ii) t) 
z o 


AI (1-cos2o 3 t) 

y o 

Al (2w t+eos2io t) 
z o o 
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For X offsets, the momentum remains cyclic and bounded with the 
same magnitude as the ideal case with the momentum values shifted 
E x 1T1 the ort:U: * For Y and Z offsets, however, a secular compo- 

nent of momentum appears which dominates the cyclic momentum for 
a 1 offset about either axis. Although the momentum peaks at 
0 ) o t- 7 r /2 and 3ir/2 for the ideal attitude, the addition, of a Y or Z 

offset of 1 ° or greater shows only local maxima or an inflection 
point moving from these positions toward oj^t^ir with the peak mo- 
mentum occurring at eo 0 t*= 2 ir with this value representing the net 
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accumulation per orbit as well. For Y offsets the peak momentum 
and momentum to be dumped per orbit is 1,200 n-m-sec per degree 
offset. Tor Z offsets the corresponding amount is 1,160 n-m-sec 
per degree offset. For offsets about both Y and Z the momentum 
accumulation is the root-sum-square of the per axis amounts. 

Figure 5 ~3 is included to indicate the effect of the princi- 
ple axis offsets on the momentum history if the vehicle geometric 
axes are used for XPOP. For no offsets, the momentum history is the 
same as the R, component of this figure. When the principal angle 

hffto^v °f h ea “ tio " S 5 " 1 * 5 ~ 2 ’ and 5 “3 are included in the momentum 
istory, the Y axis momentum dominates due to the large e . The 

momentum magnitude H is a maximum for <i> c t=2ir as described^ the 

previous paragraph. 


. Froi,: ^formation presented in reference 2, the peak aerody- 
namic torque for XPOP attitudes is on the order of 1.4 n-m with 

rnnn tranSlent: ™ mentum due tG aerodynamics being approximately 
1,UUU n-m-see. The net momentum accumulation per orbit is con- 
siderably less than this amount. 


3.2.2 X Axis in the Orbital Plane (XIOP) - The XIOP vehicle 
orientation has the advantage of allowing instrument pointing any- 
where in the celestial sphere without the need for gimballed pay- 
loads with the attendant disadvantage of more severe torque levels 
and momentum accumulation. With the vehicle X axis constrained to 
lie in the orbital plane, there always is a rotational degree of 
freedom about that axis which will be defined as 0, --ir<g<Tr the 
angle by which the vehicle Y axis is rotated out of the” orbital 
plane. . The relation between the vehicle axes and the orbit is 

ST 1 .-? f±gU ^ e 5-4 * With reference to that figure it can be seen 
that the coordinates of the unit vector along the local vertical 
resolved into vehicle coordinates is : 
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If an offset £ ~( £ x » e y »e 2 ) exists between ideal and actual vehicle 

attitude the components of the local vertical unit vector can be 
written: 
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Substituting the components a ± from equation 5-12 into equation 5-4 

then yields the gravity gradient torques for the ideal XIOP atti- 
tude. 


It. 3 2 

T = — at 
gg 2 o 


-r sin2(3 (cos2u) t-l)Al 
1 ox 

-Al sin(Bsin2tJ t 

y o 

Al cos0sin2a> t 

Z O 
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Substituting the appropriate numerical values for the differences 
of inertia and the orbital rate then gives : 


T = 
gg 


Q.l98sin2(B(eos2w t-1) 
o 

-12. 21sin$sin2u) t 
o 

11.81cos3sin2co t 
o 


(n-m) 


(5-15) 


Examination of equation 5-15 then indicates that the maximum gravity 

gradient torque occurs for sin2m t=+l, i.e., to t= ^ n=l 3 S 7 

O — 5 * O 4* 11 

further, the maximum torque about the vehicle Y axis occurs with 

| sing (=1, $=+ ~ with the value T gg =(0,+12.21,0) T . The maximum 

torque about the Z axis occurs for |cosg|=l, 8=0 ,+yr yielding at 

torque T gg =(0,0,-bll.81) T . The peak torque about the X axis occurs 

for | sin28 1=1, 8=± ~ yielding T gg =(dj3.198,+8.63,+8.35) T with 

a magnitude of 12.01 n-m. The peak gravity gradient torque is thus 
approximately 12 n-m for the ideal XIOP attitude. 
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The momentum history as a function of time can be obtained 
from the time integral of equation 5-15 as: 



89.17sin2|3(sin2co t-2to t) 

o o 

5,498sin£(cos2w t-1) 
o 

5,320cosg(l-cos2aj t) 

G 


(n-ra-see) 


For one orbit ((D^t-2-n-) the accumulated momentum is 


(5-16) 


-4-rr(89.17)sin28 

0 


-l,120sin2|B * 
0 

0 


0 

u J 


(5-17) 


This is clearly maximized when |sin2g|=l which is the value of g 
for which peak torques develop about the X axis. With the prin- 
cipa X axis in the orbital plane, the net momentum accumulation 
per orbit is approximately 1,120 n-m-sec. 

For peak momentum, the Y and Z components of equation 5-16 
dominate and both peak when cos2to Q t=-l. Since the X component 

contains the secular term 2^, the peak momentum must occur near 
w 0 t= yielding: 


-3^(89 . 17) sin2(3 


* 

-840. 4sin2£3 

-2(5,498) sing 

= 

-10, 996sin(3 

„ 2(5,320)cos8 


l0,64Ocos(3 


For peak X torque |sin2g|=l, yielding a peak momentum of 
H p =(+840.4,+7776,+7524) T with a magnitude of 10,820 n-m-sec. For 
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peak Y torque, |sin0|— 1 with H = (0, +10996, 0) 2 ' and peak Z torque 

P 

occurs for | cos 3 1 =1 with H p -(0,0,+10640) T . Thus the peak momentum 

stored in the CMG cluster is 10,996 n-ra-sec , occuring for (3=+ ~ 

for which there is no net accumulation about the X axis. A typical 
momentum history for one orbit is shown in figure 5-5 for which 
3=-ir/4. 


With any pointing offsets e, the components of equation 5-13 
must be used in the torque equation 5-4 giving the gravity gradient 
torques with an offset 



*22 
(a -a )AI 
z y x 


a a AI 
x y x 


-a a AI 
X z x 

0 -*• 2 
T =T +3 ix) £ 
gg gg ox 

-a a AI 
x y y 

+3oj 2 e 
o y 

(a 2 -a 2 )Al 
x z' y 

+3w 2 £ 
o z 

a a AI 
y z y 


a a AI 
x z z 

- ~ 


-a a AI 

L y Z Z j 


2 2 

(a -a )AI 
y x' z 


Substituting the appropriate values and adjusting the coefficients 
for a per degree offset: 


0 

T =T +e 

gg gg X 


+£ 


0 . 006 9cos 2(3 (cos 2 o) q t-1) 

0. 2130cos£3sin2(jd t 
o 

2061sin(3sin2w t 
o 


0. 0Q69cos(3sin2a) t 
o 

- . 1064 [ (l+cos23)+(3-eos23) cos2u) t] 


-.1030sin23(cos2aj t-1) 
* o 


+£ 


0. 0069sin3sin2tja t 


~.1064sin23(eos2co t-1) 
o 

- . : 03 0 [ ( 1-co s 2 3) + (3+co s 2 3 ) c o s 2 [U q t ] 


(5-20) 


>-19) 



degree, in actuality approximately 1 percent increase in neat 

per degree offset per axis. increase in peak torque 

Performing the time integral of ea nation 70 ** t, • , 

men turn history: b equation 5-20 to obtain the mo- 


0 ->- 

H „=H +e 

88 gg x 


3.112cos2(3(sin2u) t-2o> t) 

o o 

95. 95cos(3(l-cos2w t) 
o 

92 . 84sin(3(cos2(a t- 1 ) 


+£ 


3 . 112cos|3(1-cos2<j 0 t) 

o 


-47 . 96 [ Cl+cos20) 2oi t+(3-cos23) sin2cn tl 
° o 

46.41sln23(2uJ t— sin2u) t) 

o o 

! 3 . 112sim3 (l-cos2ia t) 

j ° 

+£ ; 47 . 96sin23(2to t-sin2w t) 

J o © 

-46 . 41 [ (1-cosG) 2ii) t+C3+cos23) sin2a) tl 


To determine the effect on the dump requirement. 


(5-21) 


a) t=2Tr: 
o 


| 

-+0 ->- i 

H^=H ,+ ' 
d d , 

i 

1 

L 


0 

0 


1 

r 


r 

1 

! o 


0 

e + 

X 

| —602. 8(l+cos2$) 

e + 
y 

602. 8sin23 

J 

583 . 2sin2g 


-583. 2(l-c0s2J3) 


(5-2 


For |sin23|-l, which gives the maximum momentum to be dumped for 
the ideal case, H d =(U2.5, ± 602.8,+583.2) T with a magnitude of 
1,399 n— m— sec, for either Y or 7 nf foa+n • 

momentum of 25 percent per degree offset ’about’ "theTo^Y^ef 
Considering an offset of 2.31° about Y and .085° about Z Sat‘ 
is, maintaining the vehicle X axis in the orbital Z ’ that 

to the principal X axis yields H =(1120 5 1443 7 GP f° Sed 

magnitude nf 7 o.m , d '> J ~ L ^ u ‘0 > i44j. 7,1396.8) * with a 

magnitude of 2,301 n-m-sec to be u dumped per orbit. 
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In order to determine the effect on the peak momentum state 
it is necessary to evaluate equation 5-21 for iu t=3-rr/2: 

o 



-29.33cos23 


6. 223cos3 


6, 223sin3 

-+o 

H =H + 
P P 

191. 9cos3 

E + 
X 

-452.1(l+cos3) 

e + 
y 

452. Isin23 


-185. 7sin3 


437 .4sin23 


-437 .4(l-cos3) 


E z (5-23) 


for (3 +jr, corresponding to the peak momentum of the nonoffset case: 


f0r V ®^ = (869.7 , 7911, 7655 ) T , 
E y : =(844.8,8548,7524) T , 

E z : =(844.8,7776,9270)^, 


Hp=ll,043 n-m-sec (5-24a) 
Hp = ll>4l9 n-m-sec (5-2 4b) 
H p =ll,383 n-m-sec (5-24c) 


Thus, there is an increase in peak momentum of 2.1 percent per decree 
offset about X, 5.5 percent per degree offset about Y, and 5?2 percent 
P agree offset about Z. Again considering the case of vehicle X 
axis in the orbital plane with a Y offset of 2.31° yields l Mo- 
mentum value of 12,194 n-m-sec. 7 3 peak mo 

From data contained in reference 2, the peak aerodynamic torque 

f ° r *“* With the Peak m ° meiltum bei “S 

indicatel\hf rf UVer ^ R T irem£n ; t -- " While baseline information 
the rS e rOSS Vehlcle I,laneuvers wil1 not be performed with 

the CMG system, it is evident that momentum desaturation will re- 
quire vehicle maneuvering under CMG control. The effects of this 
maneuvering can be significant in three ways, i.e., 1) the torque 
level required to accelerate and decelerate the vehicle as the man- 
e “ Ver 13 established and completed, 2) the momentum that must be 
absorbed by the CMG system during the maneuver as required by con- 
servation of angular momentum, and 3) the increase in gravity gradient 
momentum accumulation during the maneuver. gravity gradient 

sat„r*M V1 ° US Sbudies ( refe rence 2) have shown that momentum de- 
saturation can be accomplished by means of small angle maneuver 

Tn W l ln W ? 1Ch 5 he m3Xmum maneuver is no greater than 10 degrees. 

general at least two maneuvers will be required for XPOP de- 
saturation and no more thna 4 for XIOP desaturation with the sum 
of all maneuver angles in either case probably not exceeding 20 
degrees, For the 500 km orbit, the period is 5,676 seeonds^nd 
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postulating a minimum desaturation period of one-third orbit or 
1,892 seconds, no more than 20 percent of the desaturation period 
or some 378 seconds should be taken by the maneuvers. Thus a maxi- 
mum desaturation maneuver rate 0^=20/500=0.04 deg/sec seems a 

practical working value. Referring to table 3-1 it can be seen 
that the vehicle axis of largest inertia is Z^, and all following 

computations will be based on maneuvers about Z as the worst case 

in terms of torque requirements and momentum exchange. 


Conservation of momentum requires that the sum of vehicle mo- 
mentum and CMG momentum must remain constant when no external torques 
act on the system, which is nearly the case during a maneuver since 
the gravity gradient torque is quite small. During experiment point- 
ing the vehicle rates are essentially zero, thus the vehicle angular 
momentum is zero. When a maneuver occurs the CMG system must experi- 
ence a change in momentum AH=-Ju) where J is the moment of inertia 
about the axis of rotation and w is the rotation rate in order to 

satisfy the conservation law. For m =.04 deg/sec: 

z ° 


l* H z |-I zz kj =5,315 n-m-sec (5-25) 

To examine the torque requirement to perform the maneuvers an 
ideal maneuver strategy consisting of an acceleration to achieve the 
limiting maneuver rate, a coasting rotation at the maximum rate and 
a deceleration to zero rate with acceleration and deceleration times 
equal will be used. This is a time optimal strategy for equal torque 
applications on acceleration and deceleration. If the torque level 
of the CMG is defined to be T, the acceleration about the axis of 
rotation is : 

(5-26) 


The time required to accelerate to the maximum rate ui or decelerate 

from til to rest is: 

. m 


t =tu /T 

am m 


(5-27) 


The angular motion during acceleration and deceleration is: 

2 


. 2 Jw 

4 **V /2 ‘ 2T 


m 


(5-28) 


A total maneuver ip then will require an angular motion of iJ,-2Aij; 
during the constant rate period which will require a time: 
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(5-29) 


t = ^—2 Aiji _ m 

b m a) f 

m to 


The total maneuver time is then: 


, U UJ 

t =2t +t = 4- + j 

m a b to T 
m 


(5-30) 


If instantaneous acceleration is assumed, the maneuver time is simply 
i|)/a) m thus the increase in time required to perform the maneuver con- 
sidering a finite torque application is: 


At=Jm /T 
m 


(5-31) 


Using the maximum rate of 0.04 deg/sec and I then gives: 

At— 5,315/T sec (5-32) 

For an increase of maneuver time of 25 seconds (a 10 percent increase 
for a 10 degree maneuver) a torque of 213 n-m is required, while for 
20 and 10 seconds the values are 266 n-m and 532 n-m* respectively. 

A minimum torque requirement for maneuvering thus would be 200 n-m 
with a value of 510 n-m being a practical maximum as larger values 
decrease the total maneuver time only a small amount for even small 
maneuvers . 

The amount of gravity gradient torque and additional momentum 
accumulation exerted on the vehicle during a maneuver can be bounded 
by reference to the torque equation 5-19 which is a general torque 
equation in terms of the components of the local vertical vector 
S i’ 1=x,y,z ‘ Since the orbital position changes only slightly durian 
a maneuver, the a i can be assumed constant and since they are compo- 
nents of a unit vector; max (a. a j 3=0.5, max{aj-aj }=1, i*j . The maxi- 
mum torque thus arises from the e term of equation 5-19 since AI 

7 y 

Is the largest inertia difference and in this term has a coefficient 
bound of unity Thus the gravity gradient torque during a maneuver 
is thus bounded by: 


0.5AI 


T <3u £ t AI 

ggm o y y 

0.5AI 


(5-33) 
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where e y -io y t and it is a strict inequality since all three coef- 
ficients cannot attain their maximums simultaneously. The cor- 
tion°S-^ m ° mentUm b °und is found from the time integral of equa- 


t < 4 m 2 e t 2 

ggm 2 o y 


0.5A1 

2 

AI 

y 

0.5AI 


(5-34) 


F ° r £ y “m * 04 de s/see and a maneuver time of 250 seconds: 


-y 

T < 

gem 


H < 

ggm 


0.069 

4.238 

2.050 J 

8.542 

529.7 

256.3 


X <4.70 n-^m 
ggm 


(5-35) 


H ggm <588 n-m “ sec (5-36) 


- 3,3 Summary of Torque an d Momentum Requirements - The results 

of the previous sections are summarized in table 5-1, listing all 
derived torque and momentum requirements. The overwhelming torque 
rJirT 13 tha ^ f ° r desaturat±on maneuvering. A CMC cluster 

ist^r^ a minimum 0f 200 n_m Gf torc I ue about any axis 
is indicated with torque outputs up to some 540 n-m being helpful in 
reducing total maneuver times. 1 UJ - ln 

The XIOP attitude is clearly the worst with respect to momentum 
rage requirements. The peak gravity gradient, aerodynamic, mo- 

950 5 2W and 580 n ' aneWer 8ravity ^adieut requirements are 10,980, 
11° 1 5,2 ? G ’ and 580 n ~ m ~ se o, respectively, for this case. Each of 

case thS Ue ? 7 rePr r ent f/ eCt ° r magnitude and even the unlikely 
th!t ?! aka WQUld ° CGUr simultaneously, there is no chance 

that the individual vectors would be colinear. Realistically the 
total momentum storage requirement can thus be estimated by using 
the root sum square of the four values instead of the direct sum 
resulting in a minimum momentum storage requirement of 12.200 n-m- 

TtST aX±S * , IL should be noted that the largest contribution 
to the momentum requirement is that due to gravity gradient for which 
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Table 5-1. Summary of Torque and Momentum Requirements 


Peak Gravity Gradient Torque 
Increase per degree offset e 


Peak Aerodynamic Torque 
Minimum Maneuver Torque 
GG Torque - 250 Second Maneuver 

Peak Gravity Gradient Momentum 
Increase per degree offset e 


Peak Aerodynamic Momentum 
Momentum Exchange .04 deg/sec 
GG Momentum ~ 250 Second Maneuver 
Momentum Accumulation Per Orbit 
Increase per degree offset e 


0.41 

12.20 

n-m 

0.00 

0.15 

n-m 

0.15 

0.15 

n-m 

0.15 

0.15 

n-m 

1.50 

1.50 

n-m 

200.00 

200.00 

n-m 

4.75 

4.75 

n-m 

353 

10,980 

n-m-sec 

0 

68 

n-m-sec 

854 

434 

n-m— sec 

813 

407 

n-m-sec 


1.63-0 

5,290 

580 

0 

0 

1,210 

1,170 


950 

5,290 

580 

1,125 

41 


n-m-sec 

n-m-sec 

n-m-sec 

n-m-sec 

n-m-sec 

n-m-sec 


I 


270 

270 


n-m-sec 



SLe^M« d 'lf,r POne ' ,t f °* St0red ““-turn are quite predictable. 
Mae a l : Che Ca8e the lniCIal Caster momentum state can be 
retirement 6 ““ “ 8nICu<]e of the P« ks . thus easing the total 

mnma f r ° m ab ° Ve analysis it: appears that a CMG cluster havlnc a 

capabSLf or!nn Capabllit:y ° f 12,200 n “ m " sec and a Ininimujn torque 
S b T i 00 Can meet the requirements of the orbiter/pal- 

JiL Veh H Cle b ° th f ° r gravity Sradient stabilization during observa- 
tion and maneuvering for momentum desaturation. 

? S !f neralj torc l ues and momenta can be minimized for either 

the orbiKrT ^ u princ± P al 3x13 eifc her perpendicular to or in 
t' 1 p ^ ane » thus lns uring some additional system capability 
a ^ addltional disturbance events. This could be accom- 

about e the U vehic^ S± v y ^ XP ° P SinCe the ma i° r "^alignment is 
IV 116 vehicle Y axis and all experiment configurations show a 

£atetv°5 ?oT abGU 5 that thus a sim P la rotation of approx- 

imately 2.3 about the vehicle Y axis will place the X principal 

Se a rb±C r™* 1 and the ^Per^ent gimbal eo^ld compen- 

sate for that directly. In the case of XIOP, the experiments c^ 

be hardmounted with the vehicle Z axis being the pointed a^s how 

ever, a rotation about the vehicle Z axis could place the principal 
X axis near the orbital plane. principal 


5-21 


6. CMG CONFIGURATION STUDY 

The CMG cluster requirements of the previous cenH™ u j 

Momentum: 12,200 n-m-sec in the vehicle YZ plane 

Torque: . 200 n-m about any control axis 

f In addition to these basic requirements a fail operational 
cniM ^ re J Undancy re quirement is assumed to insure mission ’ 

CMG actuators presently available or in late development 

si:r S wSht UrV T d and traded by type (doutle vs s iugle girnbal) 
basJs ^ d P ° Wer requirements on an actuator and a cluster 

r ; ra li^r is 8 ^ en to both normai s er 

operation with some qualitative vibrational aspects discussed. 

* 5,1 -^g -iyPes and Type Select -in n - The two t V pes of rMP<= 
are the single gimbal CMG (SCMC) with 'one gimbal between the 1- 

two e™ h ael and the base and tile double gimbal CMG (DCMG) with 
two gimbals between the momentum wheel and the base AnictnrL 1 
representation of both CMG types is shown in figure * 6-^. ? 

Tbe SCMG has^only one rotational degree of freedom and the 
output torque is thus constrained to lie along the line 

normal to both the gimbal axis and the wheel momentum vector H 

tL iS g^bafax S i C s G th°: t:PUt tGrqUE 3X15 ’ and SlttCe * lies normal to 

gimbfrtoiqu" ^s re JbelcM ° f OUtPUt t0rqUe about the 

bv the rS!i n a out P ut torque is thus limited only 

»Lel Si"? ° f the and .omentum 

wneel Dearmgs. The gimbal torquers must only accelerate the 

SS 1 Si ”f there ale lar Sd reaction Lrques Ibont ^he 
SCMG Is^h 5 ’ re * at±vely larsd gimbal rates are possible, the 
SCMG is thus a torque multiplier, i.e., smell torques applied 

about the gimbal axis to establish a gimbal rate S can produce 
large torques about the output axis In hmr,H S 
gimbal CMG Is „ell suited £ hi^ Wrq^ applLaSoL but'L 8 . " 
a constrained momentum storage capability. » s 
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gimbals can produce an effective global rate J anywhere in the 
piane normal to the wheel momentum H w thus the^utput torque 
T D “g xH w C3n lie a «ywbere in this plane which can then be called 

optimum) momentum storage capability. iexible (i.e. , 

The momentum envelope of one CMG -in n i a , 

gimbal axis and the mveLe of l l , ? normal t0 the 

hLlLirSusLrcL^brL'ranld '* 16 indi ^ d -" S -tuator: T'ttSL 

momentum envelope, a sinele^McTf-i M SiVe 3 m ° rS ° r leSS regular 

envelop, and 

termining desaturation strategies e ^ f T problems in de- 

desaturation which tends to decrease ‘the stored 3 gradient 

could drive the cluster stored momentum magnitude 

~£.~3S 

« C to t ^ T” ±S r 0351 ” 1 " ab ° Ut «*« normal ro Jhe 

Sar a and^the ral " 1 ^ 0 l^ 3r3 °^ 3 ^ l3 ^ a ^ 3ba3 ^ b ^ la ^ 3 ^ aac ^ f ^®®® S ^ 0 ^ llr- 

The DCMG momentum envelope is snherifai -ir, t 

of m ! a single S a t d “ th limited gimbal station the* envelope 

from the envelope in both normal and f!^ ! 3y f moving awa y 

T^PMr i • nrn_i_ and failucs modes of operation 

DCMG singularity occurs only when the ind-f->Hd, " operation. 

vecrore are colineer and this situation can 

momenta 6 through^rhe°use oTl^Z °* ■*. individual 

when more tha „ 6 . Single 6 * fr “ d °" »-»*« 
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An examination of the capabilities of currently available 
actuator^o 3 t * bulated in referenee 3 indicates a maximum SCMG 
up V To50 nTs T Capabllit y of i^ 350 n-tn-sec with units 
ol dmti’v ^ ff Study * Remaining within the available 

nearly available single gimbal actuators implies the need for 

momentun^sCDrage 

,, , application (including a margin for singular:! f-v 

avoidance). The attendant mechanical complexity of at least 10 7 

Wiethe f T attendant: m ° UntlllS and ^tion problems Llong 

indLa^L thft tS'si^f le ? ^ fallUre m ° de Saturati °n avoidance 
dicates that the single gimbal CMG should not be considered f„r 

ion wmt fl ° ated f llet applicati -> all furt^r Valua- 
tion will be concerned with selection of a double gimbal actuator 

an specification of a cluster configuration utilizing DCMGs. 

Survey of Double Gimba l CMGs - Examination of the speci- 

erenceTanf or nearly available DCMG actuators from ref- 

considered Li Sr 3 Sh ° WS three devlces which will be 

considered for this application. Basic data for these actuators 

are shorn in table 6-1. The Bendix MA-2300 DCMG is the actuator 

are ir^he'lfte 1 ? 1 ™ hlcl<S ” hile the Sperry 4500 and Bendix MA-2000 
re in the late development stage. Although the MA-2300 ach.^ar 

xperienced spin bearing lubrication difficulties during the Skylab 

r h , 1Sa has . si "“ »een solved and the modLied device 
add ' r ■ f only double gimbal CMG which is space qualified In 
addition the gimbal stops on the MA-2300 have been removed 

requ^ln fot^r^n Bef ln1 ~ Th “ momentum storage 

quirement of 12,200 n-m-sec defines the low limit r>f the ire, 

Tl SperryLsOO DCMG 4 ^endiac MA-2300 or MA-2000 

ZoLl- ? 4 DCMGs are re< l ui red. In order to meet the fail 
thJV ^ C ° ncept ’ however, at least 3 CMGs are required since 
3 S tV1n COnitr01 ±S nDt available with a single actuator, thus 
bilities a ? tuators would be necessary. The total capa- 

llnstlTs trt ^ uiren,en ^ of a 11 four of the above mentioned 

, , J 81veu ln table 6-2, although the 2 CMG cluster will 

not be considered further <lu- t-o .. . SCe ^ Wl11 


opera- 


nnt . , u ai enough the 2 CMG clustc 

tional" V V f Urther due to the “Ability to provide 

tional failure mode protection. 

6-2 t'Jl — ■ Baeommandarioos - Referring again to table 

ateiy evid““ tSrin “* f" *““* «»*“*•. 1C is *— *•' 

d-iffL-J ln the norma l operating mode there is little 

LJS The in V 2 e 300 ° rrv n baS±G COntrGl -pabilities and 
ight. The MA-2300 cluster has a higher power requirement in 

tranL V ° peXation while the MA-2000 cluster shows higher 
transient requirements although both are moderate. The variable 
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Table 6-1. Candidate Double Gimbal CMG Actuators 


Ln 


L 


Stored Angular Momentum 
Peak Output Torque (n-m) 
Wheel Rotation Rate (rpm) 
Actuator Weight (kg) 

Actuator Dimensions (m) 
Simbal Freedom 
Actuator Bandwidth (Hz) 

Wheel Spinup Time (hr) 

Power Requirements (watts) 
Wheel Spin Control 
Spinup Peak 

Gimbal Torquers at Peak 
Spin Bearing Heaters (peak) 
Other (peak) 

Spin Motor Type 
Torquer Type 
Torquer /Gimbal Drive 


Status 


BENDIX 

MA-2300 

SPERRY 

4500 

3,120 

6,100 (d) 

165 

340 

9,080 

6,510 (d) 

190 (a) 

295 (c) 

lxl.Q7xl 

1.17x1.22x0.53 

A 

Unlimited 

4 to IQ 

15 

14 | 

4 


80 
170 
170 
52 

70_(b) 

AC 

Brush type DC 
Geared 


Space qualified 
Available "off- 
the-shelf" 


55 to 102 (e) 
450 
500 

Unknown 
Unknown 
Brushless DC 
Brushless DC 
Direct 


Development 


Notes: 


a) 

b) 

e) 

d) 

e) 

f) 

* 


BENDIX 

MA-2000 


1,355 to 4,070 (f) 
240 

3,950 to 11,850 (f) 
205 

1x1.07x1 
Unlimited 
5 to 10 
2.5 


50 
400 
300 
52 

Unknown 
Brushless DC 
Brushless DC 
Geared 


Development 


A1wlVr 1Ude inverter assembly weight of 23 kg, 2 required per cluster. 
Allocated to inverter assembly functions. 

Does not include external control electronics of unknown weight. 

Nominal values (4 discrete speeds possible) . 

Larger value required at peak output torque. 

Nominal values 2,700 n-m-see, 7,900 rpm. 

Unlimited gimbal freedom in modified actuator. 



Table 6-2. Candidate Actuator Clusters 


Normal Operation 

Momentum Storage Capacity (n-m-see) 

Peak Torque Output (n-ra) 

Cluster Weight (kg) 

Power Required for Spin Control (watts) 

Peak Power Required During S pinup (watts) 

Torquer Power Required at 10 ft-lb-sec (watts) (d) 
Torquer Power Required at 150 ft-lb-sec (watts) (d) 
Torquer Power Required at 400 ft-lb-sec (watts) (d) 
Single Actuator Failure 
Momentum Storage Capacity (n-m-see) 

Peak Torque Output (n-m) 

Control Law Modifications Required for Failure 


4 BERDIX 
HA-2300 

4 BERDIX 
MA-2000 

3 SPERRY 
4500 

2 SPERRY 
4500 

12,480 

660 

12,200 

960 

18,300 

1,020 

12,200 

680 


760 (b) 
320 


9,360 

495 

Minimal 


a) Run at 3,050 n-m-sec each instead of nominal 2,700 n-m-see 
b Does not include 46 kg for 2 inverter assembles. 

) Does not include weight of external control electronics 

d) Estimated from peak values. s - 

e) Includes increase in spin control power required. 

f Based on increase to 4,067 n-m-sec for remaining actuators 
g) Three axis control impossible. ’ 


820 

200 

1,600 

16 

255 

690 

12,200 

720 

Minimal 


885 (c) 
165 
1,350 
25 (e) 
330 (e) 
875 (e) 

12,200 

680 

Moderate 


590 (c) 

no 

900 
25 (e) 
330 (e) 
875 (e) 


6,100 

340 

(g) 



momentum storage feature of the MA-2000 actuator allows the cluster 
storage capacity to be returned to 12,200 n-m-sec in case of a 
single failure, however, this is not the overwhelming advantage it 
appears to be since the gravity gradient momentum peaks which 
dominate the transient momentum are highly predictable and the 
cluster momentum state can be biased to provide the necessary mo- 
mentum storage capability in the required direction. Thus, since 
the actuator capabilities and requirements are quite similar and 
the decreased momentum capability of the MA-2300 cluster in the 
failure mode is not a serious problem, the MA-2300 cluster is to 
be preferred to the MA-2000 on the basis of the space qualifica- 
tio and off-the-shelf availability of the actuators. 

Comparison of the cluster of four Bendix MA-2300 actuators and 
that of three SpeTry 4500 actuators shows that the Sperry cluster 
has higher torque and momentum storage capabilities both in the 
normal and failure modes, however, as stated above the peak momentum 
is predictable, thus the excess storage capability is not an over- 
riding consideration. The Bendix cluster shows a higher steady 
state power requirement with the Sperry showing higher transient 
power needs, although both are nominal for this application. The 
basic cluster weight is significantly less for the Bendix actuators 
and since both CMGs under consideration require additional elec- 
tronic gear, the net weight advantage for the MA-2300 cluster should 
exist even after adding the weights of the inverter and external 
electronics assemblies. Another factor to be considered in this 
comparison is the momentum wheel rotation rate of the actuators 
which defines the basic vibration frequencies arising from wheel 
imbalance. As can be seen from table 6-1, the Bendix actuator 
wheel spins at over 9,000 rpm while the Sperry spin rate is around 
b,5U0 rpm, nearly 30 percent slower. In general as the vibration 
disturbance frequency increases, the effect on pointing accuracy 
due to structural propagation of the vibration decreases. Another 
factor to be considered, however, is the fact that although the 
Sperry actuator has a lower vibration frequency, there is one less 
rotating wheel in the cluster which could tend to compensate for 
the lower rate. The effects of vibration must thus be left open 
pending some study of the GMG cluster /pallet structure interaction 
and the net effect on pointing. 

In summary, the comparison of the Bendix MA-2300 and Sperry 
4500 clusters shows a slight weight and probable better vibration 
characteristics for the Bendix cluster while the Sperry cluster 
has somewhat lower power requirements. Since detailed cost in- 
formation is not available for the Sperry actuators, this factor 
cannot be evaluated. In the absence of any clear-cut disadvantages 
for the MA 2300 actuator, this actuator is to be preferred due to 
its availability and space qualification. 
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p ,. 6 *L luste r Coilfi fi»ration - With the selection of the 

Lendix MA-2300 CMC (the modified Skylab ATM actuator) for this 
application, the next problem to be considered is the mounting 
c nfiguration of the four actuators required in the cluster. In 
general , the mounting arrangement is not critical for a double gim- 

straSts W of gimbal freedom aad ^e purely hardware con- 

° f Physical clearances, torque paths, etc., will be the 

P _ ary consideration in determining the actuator locations and 
^rientations. This aspect of the configuration is treated in de- 
ta1 ^- 111 volume III of this final report. 

For the purposes of implementing a simulation of the total 

™L\ ter ~ Pal ^ et ~ SUSPenSi011 dynamics with detailed shockmounted CMG 
” d ^’, a ‘: luster configuration was required prior to completion 
of the hardware studies. Since orientation is not a critical prob- 
em within the scope of the overall study, a representative con- 
figuration was selected for use in the simulation studies. 

. The simplest arrangement for any number of CMGs would be an 
c ° nfdgurat:Lon wifc h the output axes of all the actuators 
J- ea f lng also to a relative simple CMG control law (gimbal 
that^r 5 ^ : 1 HQWeverj Placement in this orientation would require 
in 1 h conditions on the gimbal angles would have to vary 

Sit a actuator to avoid starting with all wheel momentum vectors 
aligned (a singularity condition, actually a saturation). In order 
to avoid this problem while retaining a simple configuration, a 
rotated m-line configuration was selected in which each CMG is 
rotated by m/2, n-0,1,2,3 about the base X axis as shown in fig- 
ure 6-2 This configuration is also a worst case in the sense of 

actuator bandwidth since operation will be near zero 
gimbal angles for all CMGs. 

6.6 CMG Control law - With the selection of an actuator clus- 
ter configuration to be utilized in the simulation study of the 
loated pallet, it is necessary to derive a CMG control law which 

rir 1 1 g£nerate a seL of g im bal rate commands that will result in a 
CMC cluster torque output equal to the desired command torque. For 
the purposes of this study there is no requirement to include sin- 

CMG 3 'ei2™rf an ^ l0glC and fa±1Ure ra ° de iterations to the basic 
CMG control law, however, these items could easily be added to the 
law derived in this section. 

The CMG control law to be used is a pseudo-inverse law which 
is optimized and decoupled and must be implemented digitally. The 
eight CMG gimbal rate commands are computed by minimizing the sum 
pL t ^ ie Squares of the commands with the constraint that the ideal 
CMG output torque equals the command torque desired. The extremlza- 
tion is carried out by the technique of Lagrange multipliers. 
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The pseudo-inverse steering law will be derived under the 
assumptions that all wheel momenta are equal and that the pallet 
is inertially held, i.e., all pallet rates are assumed zero. With 
reference to the recommended cluster configuration, the CMG output 
(P') axes are located with respect to the pallet (P) axes as shown 
in figure 6—3. The transformations from the CMG output axes (P'j, 
j=l,2,3,4) to the pallet coordinates (P) are: 


fp T p'lJ 


0 0 
1 0 


0 0 


(6-1) 


[ p T p. 2 J 


10 0 
0 0-1 
0 10 


( 6 - 2 ) 


f p T pi 3 ]- 


10 0 
0-10 

0 0-1 


(6-3) 




1 0 

0 0 

0 -1 


0 

1 

0 


(6-4) 


For each actuator the transformations from output (P'j) to outer gim- 
bal (Cj) and outer gimbal to inner gimbal (Aj ) coordinates are: 


[ cj T p , j ]: 


eos5 sin6 . 0 

3j 3j 


-sin6 , 
0 


eosS 


3j 


(6-5) 


F T 1 = 
L Aj i Cj J 


0 

COS if) 


Ij 


sintS. , 
ij 


-sinS . cosd n . 
Ij lj 


( 6 - 6 ) 
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where the gimbal angle 6 refers to the ith gimbal (i=l inner, i=3 
outer) of the jth CMG (j=l,2,3,4). 


With the assumption of small vehicle rates and ideal actuators, 
the total angular 4 rate applied to the wheel of the jth CMG is due to 
the gimbal rates 6^ and 6^ , explicitly: 



• < 


■ « 

0 


• m 

• 

6 lj 
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0 

cr 

V ln4 l 3 

0 

* * 


6 

. 3 '1 
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In the inner gimbal space the wheel momentum H A j=(0,h,0) T and the 
external torque applied to the wheel is the inertial derivative of 


DH 


— 6- .cosS 1 . 
3j Ij 


7rr J ~ -H ,-ko .xH. .=11 
Dt Aj Aj Aj 




The individual CMG control torque is simply the negative of the 
external torque acting on the wheel. Transforming to the CMG 
output space: 


-ntp'jVtcjV 


.cosiS- . 
3j lj 


-^IjSiniS^jSinS^j+S^eosSyCosd^ 

0 

• 

=H 

5 lj sinS ij cos6 3j + Sj COS!S 1: j s incS 3 . 

. - 4 5 13 

” Mi 


-V osa ij 


The total control torque exerted on the pallet by the CMG cluster re- 
solved into pallet coordinates is then: 
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or: 


4 

T CP = 2 !!! ^P T P ? j^ T p*j 

i=l 




(6-9) 


(6-10) 
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where: 


( 6 - 11 ) 


¥< S ii>«iA 3 .V T 

V<*31-*32-*33.*34> r 


(6-12) 


sinfi 11 sin6 3l sinS^sinfi^ -sinfi^sinfi^ -sin-fi^sinfi^' 

ll COSl5 31 cgs6 12 “ sin<S 13 COS(S 33 ~cosfi 14 (6-13) 


-costS, 


siru5 12 cos<S 32 


-sinS., cos<5„ . 
14 34 


cos6 n cos6 3i eos6 12 cos6 32 cosS^eosd^ oosS^cosS^ 


IB]- co S 6 llS in6 31 


-cos6 13 sln6 33 


cosfi 12 sin<S 32 


-cos<S sin6_. 
14 34 


(6-14) 


In order to derive the. control law it is necessary to find the gimbal 
rate commands 5 1{ , and 6^ such that the control torque T is equal 
to the command torque T’ while minimizing the generated command rates. 

lnd£X P *** iS - 


P 2 5] (6 liC +6 3iC )= I ( ^lC*\e + ^3C*^3C } 
-,* — 1 


(6-15) 


constraint * is ^simply^ ^ ^ pr ° dUCt {lnner P^duct) . The 

V*C =G (6-16) 

The optimization is carried out by using the Lagrange multiplier 

A (A 1’ A 2* A 3 ) ' the Lagrangian adjoint equation using equa- 

tions 6-14, 6-15, and 6-16: 


L=P-t.(f cp -*c> 


L " hKc'Kc + Kc'Kc ) - t ' { lU A ]iic+ H [ B ]i 3c -T» } 


(6-17) 
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To minimize L, the partial derivatives 


\ 3 l 

— , and — - are set 


3 *1C 3 *3C 


9A 


to zero, resulting in 11 equations in the 11 unknowns (45 4o 

liC 3±C* 

and 3A ± ). The notation of the derivative of the scalar function L 

with respect to a vector is understood to be the vector whose com- 
ponents are the derivatives of L with respect to the vector components, 
i.e., for any a, & £(|^-, iiu .. 

sj 8 “l S “2 3 V • 

The derivatives are reasonably straightforward perhaps with the 

* 

exception of the form ~—(Z- [C]t ) where [C] = [A] for i-1 and [B] 


9^ 


for i“3. Expanding: 


1C 


• • 

[C ^1C =A 1 (C ll^ilC +G 12^i2C +C 13^i3G +C 14^i4e ) 
+A 2 (C 21^ilG +C 22^i2C +G 23 5 i3C +C 24 S i4C ) 
+A 3 (C 31^ilC +C 32^i2C +C 33\3G +C 34^i4C ) 


C ll A l' fC 21 A 2 +€ 31 A 3 


, CA- CC3^ iG ) ! C 12 A 1 +C 22 A 2 +C ^? A 




iC 


12 1 22 2 32 3 
C 13 A 1 +C 23 A 2 +G 33 A 3 
C 14 A 1 +C 24 A 2 +C 34 A 3 


J 


where the superscript T indicates the transposed matrix. Thus 


[A] T A=G 

3 *ic 


(6-18) 


~~ =$ 3c -H[rf%0 

3l 3G 


(6-19) 


% =T c" H 

9A L 


1C 


3C 


=0 
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Substituting 6-18 and 6-19 into 6-20: 

T^H[A]lI[A] T t+H[B]E[B] T t 
=H 2 {[A] [A] T +[B][B] T }t 
Defining: [G]»[A] [A] T +[B] [B] T 


(6-21) 


( 6 - 22 ) 


fG] is a 3x3 matrix and is symmetric by examination of equation 6-22, 
further, since it is square, equation 6-21 can be solved for t: 


V G] " 1? c 

H L 


(6-23) 


Finally, substituting into equations 6-18 and 6-19 and rearranging: 


V itA] T t Gr^ 


(6-24) 


(6-25) 


mu summary , the CMG control law must be implemented digitally. 

he gimbal angle state is required to form the matrices [A] and [B] 
as defined in equations 6-13 and 6-14, respectively. The matrix [G] 

is formed using equation 6-22. Then, given the command torque T\ 

the gimbal rate commands are computed from equations 6-24 and 6-25. 

It should be noted that the basic steering law presented in 

this paragraph can easily be altered for changes in CMG mounting 

arrangements. Once a cluster configuration has been defined, the 

output to pallet transformations can be determined and these will 

used"?; ! qu % tiGns fi 6 : 1 thrGu ^ • These new transformations are 

TBl aS generate new forms of the ^trices [A] and 

|BJ and the balance of the computation remains as given above. 
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7. 


CONTROL MOMENT GYRO (CMG) MODEL 


(ATM) T CMG CMG m! ed - 11 th±S St " dy ±S Skylab A P°Uo Telescope Mount 

v"us studief^M? l,Ue “£“5 ° f thiS Study ’ * h “ to pre- 

, studies utilizing the ATM CMG model is that rh P run -r ^ 

vious d s t°i the fl ° ated pallet throu e’i shockmount isolators. i/pre- 
, ri i adles associal: ed with the ATM CMG, the CMG is considered to 
ratei 8 ^ 17 !"° Unted t0 the vehicle * Also in previous studies vehicle 

therefore^M^base^ates^er^assumed to^TzIro" ^ ^ ”*"• 

sprinps nC ?he n ha^ e SCUdy the base of the CMG is mounted on 

P gs, the base rates have values comparable to the gimbal rates. 

The primary purpose of developing a CMG model for this study 
is to determine how the CMG base rates (^) , caused by shockmount 
dynamics, change the six-mass CMG dynamic model. 

to be 7 rig l f °with a th e - S - ~ A11 . structures the CMG are considered 
co De rigid with the exception of the vea-r f-t-n-ir. • 

the compliance between the loner gimbal and the CMG retort *” 

7 -2 Sh ockmount - The CMG la » 0l ,nted on shockmounts along each 
axis. The purpose of this shockmount system Is to Isolate thf 
frequency vibrations within the CMG from reachinc the n n t " 
Shockmounts have a bandwidth of 10 H f “SffrSofi 

wfdth „f ffcMG^Ihf °“,f e CU ° SyS “” " U1 bA b ° l™«° fe baL 

be fremipnrv MG ‘ Thls j wi11 rec i uire that the CMG /shockmount svf lT1 

«/shS^“ CO lncreas ‘ the band ” ldtb bba overall ' 

switch’h 2&- ™***** .- The CMG friction is modeled classically to 
switch between static (stiction) and running (coulomb) faction! 

Ideally, this switch occurs whenever the motor rate 6 . (i=l,3) or 

gimbal rate <S ± is driven to or away from a zero value. The following 

gimba! torque motor friction , ^utej gimbS 8 tfrque^ moto^ f 

V“~ ^ teXISfJf 

stiction and equal to the running level T M* durinfcoulomb 
friction. SJ - g 
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During stiction, the friction torque T -T , achieves what- 
ever instantaneous value is required to maintain the inner gimbal 
fixed relative to the outer gimbal, i.e., ^=0, so long as this 

torque never exceeds the maximum static level S . Thus, during 
stiction ® 

0) -Ui 
ax cx 


and 


U) =0) 
ax cx 


where 


“ax inertial inner gimbal angular rate about inner gimbal 
x axis 

“cx = inertial outer gimbal angular rate about inner gimbal 
x axis 

The torque required to maintain this condition is therefore obtained 
as follows: 

-T t=T ,-T -J u 
sgl gl wx all cx 


where 


T g l ~ in-ner gimbal gear train output torque 

T wx = wheel momentum reaction torque about inner gimbal x axis 

J all = inertia of inner gimbal and motor rotor assembly about 
inner gimbal x axis 

As a practical matter, each CMG friction model was implemented 
on the analog portion of the hybrid computer as shown in figure 10-2. 
e friction torque switches from a running condition to a 

static condition whenever 


and 




<S 

S 
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where is a small value slightly larger than the noise level of the 
computer. Once in a static condition, the friction torque T f ^ switch- 
es to a running level only when ° 


This insures that only the required breakaway torque will initiate a 

running condition and not 5 noise. The value for T is monitored 

x sgl 

and tested continuously' in order to implement these switching deci- 
sions . 


Rather than feed T fgl -T sgl to the inner gimbal inertia integrator 

during stiction, and thus introduce the possibility of drift between 
inertial inner and outer gimbal angular rates, the inner gimbal in- 
ertial rate (0^ is forced to equal (l.e., ^=0) by disabling the 

inner gimbal integrator and causing the initial condition to equal 
^cx at any instant of time.- Once switching to the running condi- 
tion has occurred, the integrator is then freed to follow the input 


T t -T , =T -T +T 

gl wx fgl gl wx Rgl 


Furthermore to avoid gear train stiffness integrator drift, the gear 

tram stiffness integrator is placed in hold when both motor and eim- 
bal are stuck. 


• f ' iaally » t0 avQ id an algebraic loop in implementing the outer 

gimbal static friction, the following calculation is made. The 
expression for T „ is 
sg3 


=-T _+c6-T +J 


sg3 J 'g3 ' '" u l' L wz ' J c33 U bz 


where 

T sg 3 “ outer gimbal static friction torque 

^g 3 _ outer gimbal gear train output torque 

T wz = wheel Ii,omentum react ion torque about inner gimbal z 
axis 

J c33 = ec l ui valent outer gimbal inertia about outer gimbal z 
axis 
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angular rate about outer gimbal 


01 , 


bz 1=3 inertial CMG base 
z axis 


Now, w bz can be evaluated from the z axis base equation of motion as 


follows : 


where: 


w bz = ~ 7 — “(T ,,+T ,+T ) 

b33 F ’ 3 s S 3 P-z 


J b33 = CMG base ^ertia about base z axis 
Vz = torque exerted on pallet about CMG base z axis 

Substituting into the expression for and solving for gives 

J. 


T =-T - — 

s g 3 g3 J. 


J c33 


4-T f 

b33 c33 


T + — 
pz J. 


b33 


+T 1 

b33™c33 


cos<S,T 

1 wz 


This expression for T sg3 is used in she friction logic on a continu- 
ous basis, and is also substituted for T in the shockmount dynamics 
during outer gimbal stiction. 

space) 1 ”)™ 6 Sl ’ Sckm ““" ted ”«“ ;Th irt” e „Sel c“rdinate 

“r “ r srei^’f^iBis^r 

end 5) the pallet end of the shockmount (P'-speca) Fieure iTJ’ 

CMG ££* , 0f these -tdiLte ; ? sXr ^„tt 

gimbal angles Cq.c,) are in the zero position and the two spring 

equivalent systems (Y ± ,e si ) are in the relaxed state, than all five 
systems are aligned. 

transformations defined? “ ““ «- -^opriete 

right-hand system is completed by the z axis. 4 U 
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H-in*.- 7 ’ 4 ' 2 - -^ ner timbal Coordin ate System fA-sp arel - This coor „ 

system i^def - 1S H flXe ^ ±n ^ inner simbal frame of the CMG, The 
bal tn efiaed Wlth the x axis along the axis of the inner eim- 

Z axis ^r eCt ° r ° f ^ r ° t0r the 

Compliance springs (y.) provide the connection between the 
th !°! the inner eimbal. These springs provide support along 

— sir = ;; l£*- 

the 2 axxs is defined by y y Deflections are considered to be 

small angle and the transformation from A-space to W-space is 
given by equation 7-1. space xs 


y 3 0 


[ I ]= -Y 1 
w a '3 ± 


0 -Y 1 1 


10 0 


Y x « I 0 1 0 


(7-1) 


0 0 1 


is f ZL a' 3 Gimbal Coordina t e System (C-space) - This system 

fixed in the outer gimbal defined with the x^axis through tv,* 
mer 8 ^t. r1u „, the 2 ^ through che 

and the y axxs completing the right-hand coordinate system The * 

STJ.S5E 'fir ce) “ d °T r slmbal <►*-> « rS:J 'S 

the rotation of the inner gimbal through the angle 6 Equation 
7-2 gives the transformation from C-space to A-space. 


t a T e J“ 0 cos^ sin<5 1 

0 -sin5^ cosfi^ 


(7-2) 


, 7 : 4,4 . Base C oordinate System (g-snace) - The base system 
fixed m the CMG base with the z axis along the centerline of H 

When 5 3 0, the base (B) and outer gimbal (C) systems are aligned! 
The. transformation relating the B and C system is given by equation 
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I 1\]= 

C D 


cosd. 


-sind. 


sind. 


Cosd 


3 v 
0 1 
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7 ‘ 4 ; 5 Sh p . ckmourit Coordinate S y stem (P’ -space) - The shockmount 
coordinate system is fixed on the pallet end of the shockmount. This 
system is related to the pallet system (P-space) by the transforma- 

^ u in SEGtion 6 * 6 ‘ A1sq this system is essentially aligned 
Pn *\ the ^, se ^stem. displaced only by the small angle 0 rotations. 
Equation 7 4 gives the transformation from the P'-space to B-space. 
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, 7 ' 5 . dynamics - Dynamic equations are written for the rotor 

ty-dynamics) , inner gimbal, outer gimbal, and CMG base. These equa- 
tions are used to develop the CMG rate loops and the torque output 
exerted on the pallet from the CMG. All product of rate ‘terms are 
considered to be small and are not used throughout the development. 


7 ‘ 5,1 Itiner Ginibal/Roroc Complia nce (y) Dynamics - In this sec- 
tion the equations are developed which relate the rotor (W-space) 
dynamics to the inner gimbal. The angular momentum of the rotor 
(li ) written in W— space is 
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where and are the moments of inertia of the rotor about a dia- 
meter and the spin axis, respectively, w w is the inertial rate of the 
wheel frame and w, is the spin rate of the rotor. Torque (T ) 
exerted onto the wheel frame written in W-space is 
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+ dH 

T = ~ -ho xii 
w dt w r 


•" — 


• 


— 

T 

wx 


J ,w 
d wx 


—to H 
wz r 

T 

wy 

= 

J (to +U) ) 
r wy s 

+ 

0 

T 

wz 


• 

J.tU 
d wz 


u H 
wx r 
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ti- c P f dUGt ° f rate 5 erras have been neglected because torques from 
this source are small compared to the gyroscopic torques. 


Torques exerted onto the wheel frame (i.e., T i= x ,y jZ ) are 
from the y compliance springs and are 





T 

wx 


~ lC yl T l"V^l 

T 

wy 

ss 

“ K y2 Y 2^ e y2 y 2 

T 

wz 


“ K y3 Y 3~ 6 y3 Y 3 


where Y^ is the displacement between the wheel frame and the inner 
gimbal frame, and 


• 


— « 

Y 1 


(0 — 0) 
wx ax 

• 



Y 2 


& 

3 

1 

a* 

• 



Y 3 

L 


to —to 
wz az 


where u)^ is the inertial rate of the inner gimbal. From equations 
7-6 and 7-7, the following relationship exists 


J d“wz = ' K Y 3 Y 3" e r3V"r‘wx 


Using these relationships and equation 7-8, 
dynamics in block diagram form. 


figure 7-2 gives the y 
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7-5 ' 2 - nner Gimb al Dynamics - Angular momentum (H ) of the 
inner gimbal written in A-space is 3 
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J a22 
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J a22“ay 
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where J g . are the principal axis inertial terms for the inner gimbal. 
The inner gimbal dynamic equation written in A-space is 


T 

ext 


It ^ a 


where is the sum of all external torques acting on the inner 

gimbal. External torques acting on the inner gimbal consist of 
two components, the torque exerted by the outer gimbal (T ) onto 

the inner cymbal exerted frora clie wheel fJame onto 

the inner gimbal. This second torque will be the negative of the 

torque T . Therefore 
w 


f T ^ 

ax 


T 

wx 


r T * 1 

J . - (0 
all ax 
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ay 
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wy 
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J a22^ay 

T 

az 


T 

wz 

: 

J a33 W az 
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with the product of rate terms being neglected. 

7 ' 5 / 3 r— C Gimbal Dynamics - Equations representing the outer 
gimbal dynamics are developed in a similar manner as thosf for the 


inner gimbal. External torques acting on the outer gimbal are T 
(torques from the base onto the outer gimbal) and the reaction 
torques from the inner gimbal -T & transformed to C-space. 


The angular momentum ($.) of the outer gimbal written in C~ 
space is 
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(7-11) 
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Wher ^ J ci are the P rinc ipal axis inertia terms for the outer gimbal 
and oj g is the inertial rate of the outer gimbal. 


. A eain, neglecting product of rate terms and 
grrabal dynamic equation in C-space 


writing the outer 
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T -[ T ]T = 
e c a J a dt 
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Now, writing the inner gimbal inertial rate in terms of outer 

gimbal inertia rate (w.) plus the inner gimba/rate 
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in C-space. 
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Equation 7—15 reduces to 


,.10 

ell cx 


c22 cy 


•1 

c33 cz 


“ — 


T 

cx 


T 

cy 

- 

T 

cz 



1 0 
0 cS 
0 s6 


l - s6 il 


cdL 


wx 

0 


wz 


J cll 


0 

0 


T * 

c22 


0 

0 


J c33 



“• -1 


[— " -I 


w 

cx 

• 


all 1 


0> 

cy 

• 

U) 

cz 

+ 

0 

0 


— — 


— — __ 


(7-16) 


where (torque about rotor spin axis) and (J a99 -J^ qq ) 
small and 


a 22~^a33^ are CGns idered 


J cll J cll +J all 


1 

2 1 

1. 


J c22 J c22 + 2 (J a22 +J a33 ) 


J c33 J c33 + 2 (J a22 +J a33 ) 
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7 * 5 ' 4 CMC Base Dynamics - Angular momentum (H^) of the base 
written in base space is 


H b = 


J bll 

0 

0 


J biiV 

0 

J b22 

0 

V 

J b22 W by 

0 

0 

J b33 


J b33 W bz 


(7-17) 


where ^ is the inertial rate of the base frame, and J . are the 
• • bii 

principal axis inertia terms of the base. 

Writing the base dynamic equation in B-space 

_> - dH, 

Vth T = ~?T+co u xH 
b b c c dt b b 


where T b is the torque exerted on the base by the shockmount springs 

and “tb T e^ T c iS the reaction torque from the outer gimbal onto the 
base expressed in B-space. This equation readily reduces to 


T bx 
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T by 
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c« 3 0 

T bz 
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0 1 
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cx 
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b22 by 

g3 fg3 


J b33“bz 
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and the outer gimbal acceleration is related to base acceleration by 
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when product of rate terms are neglected. 


(7-19) 


7-13 



(bi ) 7 : 5 '^:l Mse Rate About the z-Axis - Base rate about the z-axis 
( bz^ ls obtain ed directly from equation 7-18 since all torques are 

readily available. 

T bz is the torc l ue exerted on the base by the shockmount spring 
and T g3 +T fg3 is the resultant torque generated by the outer gimbal 

aW e the h“ Ce b ° th ° f 5 hCse t0rc = ues are known and act directly 
along the base z axis, the following relationship is written 


T bz ( T g 3 +T fg 3 ) J b33 w bz 


and w bz and are obtained directly as shown in figure 7-3. 

x (to 7g g V~ g - Ab ° Ut th \ X * nd 7 Axes - Base rates about 

x tOJ bx J and y (to fay ) axes is not obtained as easily as was to . T 

and T cy ln ec f ua ti°n 7-18 must be obtained by combining equations 
7-16 through 7-19 as follows 
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Expanding equation 7-20 gives 
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where 
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C6 3 T wx +s(5 1 s5 7 T 
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c6 T 
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U bz J 


J *U C6 3 S 1 

J all StS 3 fi l 

d c33^3 
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figure 7 3. Z Axis Dynamics of the CMC Base 
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[J']= 


and 


T ii 
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J 13 
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J 33 


J ix- J bu + l (j ;u +j ; 2 2 )+ i( j ;u- j ;22 )c2S 3 
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J 22 =J b22 + I (J g11 +J c 22 ) " 2 (J cll" J c22 )c26 3 
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Considering only the x and y component of equation 7-21, and re- 
arranging 
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For the purposes of this study, (therefore J* ) 

negligible when compared to J' and J' ; J» and J' 

11 22’ 11 ^ uu j 22 

imated as 


V 5 A 


J all s4 3 S l 


is considered 
can be approx- 


(7-22) 


J lr J bll + 2 (J cll +J c22 +J all )+ i (J a22 +J a33 ) 
J 22 =:J b22 + l (J cll +J c22 +J all )+ i (J a 22 +J a 33 ) 
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Therefore equation 7-22 can be solved directly for co and co 
tion T 7-22 4 S1VeS the block diagram showing the solution of equa- 

7 ' 5 ’ 5 Shockm ount Dynamics - Shockmount spring stiffness K 

and damping D g _^ have been chosen to give a natural frequency of 

approximately 20 Hz with a damping of 0.15. Since the rotation 

the r "f Sy f tem 15 Sma11 an8le ’ the t°rque exerted on 

the pallet via the shockmount is given as 


-y 

* 


T P’i =K Si 0 i +D Si 0 i 


(7-23) 


where 0^-Wp and ^ and S p are the inertial rate of the CMG base 
and pallet, respectively. 

The torque exerted on the base (T^) via the shockmount is 

the negative of equation 7-23, therefore both torques are generated 
as shown an fi 8 ure 7-5. Also. since the pallet rate is 

pared no the base rate, 0 *oj 

b 

. . 7,6 £Mg/Sho ckmount Dynamic Model - The entire CMG/dynamic 

10°(fiEure 10 “ almulabed for this study is shown in section 

0 (figure 10-2) . The tachometer loop and the motor gear train 
model is shown in figure 7-6 as established in previous studies 
(references 4 and 5). Some simplifications have been made to con- 
for V D, r° mPUt:e % SP K - e ’ Wver » these simplifications have no effect 

includSe voltL 1 v F±gUrG 7 " 6 shows the tachometer loop 

includmg voltage limiters, current limiters, and torque limiters 

which were not included in the simulation model. Signal levels in 

this study were below those of these limiters, therefore these 

limiters along with the back emf gain (K^) were eliminated from 

the simulation model (figure 10-1) . 

10 2 Effects of the CMG being shockmounted are evident in figure 
10--2. Outer gimbal rates are related to base rates fay equation 
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Figure 7-6. CMC Rate Loop 
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8. CMG-PALLET CONTROL LOOP STUDIES 

The objective of the control loop studies is the specifica- 
tion of analytically determined rate, position, and position in- 
tegral gams m the vehicle control law such that the pallet control 
loop m stable has the required 2 Ha bandwidth, and such that t^e 
pallet can meet the pointing accuracy requirements. In addition 
a compensation network analysis is performed to insure that the 
a ove requirements can be met when the specified CMGs are placed 
in the control loop. 

2 !Sw the ana J- y ® eS ° f this action are based on a nominal 
T * ?^ et CGrit rol loop, the definition of gains for other 
bandwidths parallels these derivations. Indeed, gains were 
derived for several bandwidths above and below the nominal for 

ase m the simulation study which is discussed in section 10 
or this volume. 


8*^ Vgbicle Control Law - The function of the vehicle con- 
trol law is to compute a CMG torque command T q from sensed body 

motion that will compensate for unwanted vehicle motion caused by 
disturbance torques acting on the vehicle. The vehicle control 

vehicle a ^ USed tG generat ! torque CQrnman ds required to implement 
yehicle maneuvers or attitude offsets. The specified vehicle con- 

law 3 Staadard rate plus Position plus position integral 

law which in complete form would be: 

(0 B -e + h+ [Kpi ] {/^4+hdt} 

where Che gale matrices fK R ], [K p ], and [Kp p ] are diagonal, the 

vectors and 0 Q are the desired vehicle rates and positions 

which would be obtained from maneuver control logic, u> and 0 are 

the sensed vehicle rates and positions and e represents any small 
angle attitude offset command, e.g., a CMG desaturation maneuver. 

Since the scope of this study includes only pointing perfor- 
mance, all maneuver considerations will be ignored, thus in par- 
ticular <u 0 and £ will be defined to be zero. During any period of 

H 0intin8 the coatro1 3X63 wi H be held in some constant 
relationship to a set of references or ideal pointing axes. For 
the purposes of this study this relationship can be assumed an 

identity, thus can also be set to zero giving: 
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V- [ V“" [k p^- Ik p^/ 

this llZl iy kfT P ^ rfect sensors are assumed for the purposes of 

- »>p^prj£° :r^r:^ r “L“L s 

axis, thus identifying u as $ p , the control law becomes: 

v- [ vVtv/ Vt- [Kpi ]// i pdt 

In terms of the Laplace operator s this becomes on a per axis basis: 


T ci K Ri 0 Pi~ 


jr !p± „ 9 Pi . 

pi s _K pii 2 l=x »y> z 
s 


or. 


0, (B)4 ~K_ ii . fa 


0 , 


“ k ri- 


i=x,y,z 


(8-1) 


Pi 


it i. 8 ;Lu St^^ SaS ^ L °° P ' F ° r «“ !"**“« of this study 
t is assumed that the pallet is an independent body under the con 

trol of the CMG cluster, i.e., the shuttle dynamics are ignored 

they affect the pallet as disturbance torques aJpSd 
through the suspension system. Since the pallet products of i-nwp.n*. 
are quite small as compared with the moments of inertia about the & 
control axes they are neglected making possible linear single aylc 
control loop studies. Tinder these assumptions a addS^al 

assumption that the CMG cluster shows minimal cross-axis torjueout- 

ure’g-l 6 where lGGP f ° r eaCh Pallet axis ±S as shoW11 in fig- 

law (equation 8 n ? e?a '5* Ve sign GG nurring in the vehicle control 
K 8 ~P haS beetl Garr±ed around the loop to the summing 
junction where the disturbance torque T D is applied. The subscript 

i representing the axis under consideration will not be used in the 
uccee ng paragraphs for simplicity in notation. The disturbance 
torque T 0 acting on the pallet arises from torques applied by the 

suspension system. Gravity gradient torques are sinusoidal with a 

see^and^ill^h 00 bhe ° rbital perlod > or approximately 0.0022 rad/ 
alt4rS?nn f, t J ansmitted through the suspension with no 

sec, indS ^IrSin^lfs : ^ ***' 
torque with the magnitude taken as the peak Itlue oTlT^ZT^ 
the vehicle y and z axes as determined in section 5. On the other 
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a ^ d ’,i? r ? UeS dUG to Crew motiori are of relatively high frequency 
and will be filtered both by the shuttle dynamics and the suspension 
c araeteristics before acting on the pallet. The vehicle control 
law G^s) is the rate plus position plus position integral law de- 


scribed in section 8.1. The CMG torque transfer function is de- 
ined as C(s) and will be discussed in section 8.4. The compensa- 
tion G^Cs) is discussed in the following sections, 8.5 and 8.6, 


Analysis With Ideal CMG - In order to make an initial 
determination of the loop gains it will be assumed that the CMG 
transfer C(s)=l. The compensation G c (s) is thus also unity and 

the open loop transfer from figure 8-1 is thus: 


T 

G(s) = 

E 


Js 


The design approach to be followed is to determine the rate and posi- 
ion gains and K^, i.e., setting 1^=0, with the constraints of a 

2 Hz closed loop bandwidth and .707 damping on the dominant roots, 
then to add the position integral gain determining its value 

such chat the bandwidth and damping are not significantly affected. 

. 8 * 3,1 determination of Rate and Position Gains - Setting the 
position in.egral gain to zero, the open loop transfer is: 


s 


Defining : 

R *V J 

(8-2) 



(8-3) 

G(s) then becomes 




G(s)=R ^±f 
s 

(8-4) 

The characteristic 

equation of the loop is: 



A(s)=1+G(s)=1+R 

(8-5) 


s 
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ll° m equat±on 8 " 5 the system is stable for all 

ure 8 2 7h' - JZ°*a PlGt ° f equation 8-5 is given in fig* 

re 8 2, showing the double pole at the origin and the zero at -a 

For values of R<4a, the locus is a circle of radius a centered at* 

the lift .r?h^o a ^ r '^i? C “ 15 ° n nesatl '' e real “*• “ 


root is r direrr? anip j n8 C ~l 707 ’ the angle 9 is 45 ° and the desired 
r v It directly above the zero at -a for which s=-a+ja. The gain 

beiL f P ^ nt CaU be dete ™ iaed in several ways, the most direct 
being from the magnitude requirement: 


R |s+a 

R 1 2 “-t, s=-a+j a 


or > R=2a 

The closed loop transfer of the simplified 

T(s) _ RCs+a) 
T D (S) s 2 +Rs+Ra 


loop is: 


(8-6) 


(8-7) 


By definition, the bandwidth w B is that frequency at which the closed 
loop log magnitude passes through -3 db or 


i TCioi) 1 VT 

l : D (ja))l“ 2 (8-8) 

From equation 8-7: 


i (.1 to) _ R(a+jcj) 

T_ (jto) 2 

D J ' Ra-oj +jRw 


Using equation 8-8: 


R 2 (a 2 -to 2 ) T 
(Ra-to 2 ) 2 +R 2 0) 2 ^ 

D B 


(8-9) 
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R=2a 



Figure 8-2. Root Loeus of the Ideal Loop With 11^=0 


f^ ranSl ?? eq “ a 5 ion 8-9 and substituting the damping requirement 
from equation 8-6 then gives: B 4 c 


(R/w b ) 4 +8(R/w b ) 2 -4=0 

or ’ (R/^ a ) 2 =-4+2 VT 

Since R and 0J fi are both real, the upper sign must be used and for 
any desired bandwidth with a damping ratio of .707: 


R=(2 V5-4) 1/2 o)„ 

B 

(8-10) 

a=R/2 

(8-11) 

requirement ia=47r and: 

B 


R=8. 6346 

(8-12) 

a=4 . 3173 

(8-13) 


Substituting these values into equation 8-7 then gives: 
Us) _ l+s/4.3173 


V s) 


1+ VTs/6.1056+(s/6.1056) 2 


(8-14) 


Equation 8-!4 explicitly shows the dominant root damping of V?/2=.707. 
the closed loop transfer and further that the transfer includes a 

a3 ~ 4 * 3173 radians and the quadratic lag at 6.1056 radians 
The bandwidth by derivation is ay=47r=12.566 radians. 

.®* 3,2 Determination of P o sition Integral Gain - To determine the 
position integral gain K pi we now return to the complete vehicle con- 
trol law with: 


Kp 

G 1 (s)=K r + — + PI 
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The open loop transfer is thus: 


S(s) , h 5 

J 3 

s 


GCs)=R - s - + as+b 


(8-15) 


where R and a are as defined in equations 8-2 and 8-3 and: 


b% PI /K R 

The characteristic "iquation becomes: 


(8-16) 


A(s)= 1+R - s - +as+b 


(8-17) 


which for b-0 reduces to equation 8-5. Forming the Routh array for 
equation 8-17 : 


1 Ra 

R Rb 

Ra-b 0 

Rb 0 


The Routh-Hurwitz stability criterion requires that no sign changes 
of C b r a re- the firSt C ° lumn ° f the arra y» thus th e allowable values 


or. 


0<b<Ra 

0<b<37.278 


The closed loop transfer function is now: 

T( s ) s 2 -fas+b 

T (s) K 3 2 — 

R s +Rs +Ras+Rb 


- =r (b-ui )+1 am 

R(b-U) 2 )+j(j(Ra-tu 2 ) 


C8-18) 

(8-19) 

( 8 - 20 ) 

( 8 - 21 ) 
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In order to evaluate the effect of adding the position integral gain, 
the magnitude of the transfer function 8-21 is evaluated for various 
values of b at the desired bandwidth frequency of 2 Hz (4ir radians) 
using the values of R and a previously derived. 

It is evident from examination of table 8-1 that the addition 
of the position integral gain has only a minimal effect on the log- 
arithmic gain of the closed loop transfer, i.e., the bandwidth w 
, B 

is a very weak function of b. Indeed, for a selection of b=2, the 
gain is down only ,005 db from the desired -3.01 db at w=4tt. To 
determine the effect on the dominant root damping the values of R 
and 3 from equations 8-12 and 8-13 along with the values of b from 
table 8-1 are substituted into equation 8-20. The resulting cubic 
in che denominator can be factored as can the quadratic numerator 
yielding the critical frequencies of the transfer function shown 
in table 8-2. 

Examination of table 8-2 shows that the effect of adding the 
position integral gain is twofold: 1) a lag-lead factor pair is 
superimposed on the basic rate plus position control law, and 2) 
the lead and quadratic lag critical frequencies are decreased with 
increasing K^. The dominant root damping is little affected by 

the position Integral gain, however, as b increases the location; - ; 
of the position integral lag— lead pole zero pair diverge and the 
resonant peak resulting from the position gain lead and the posi- 
tion integral lead increases in magnitude. The selection of b and 
thus K pi is necessarily qualitative, however, a value of 2 seems a 

good compromise showing almost no alteration in either bandwidth 
or dominant root damping and virtually no change in the closed loop 
response of the system as can be seen in figure 8-3 which shows the 
closed loop Bode plot for b=0 and b=2. Utilizing the pallet inertia 
data from section 3, the preliminary design values of R, a, and b, 
and equations 8—12, 8-13, and 8-16, the numerical values of K^, 

Kp, and for each vehicle axis can be determined for the 2 Hz 

loop. The development for other bandwidths is similar and is not 
presented here, Gain values for all bandwidths considered are 
given in table 10-3. 


CMG Characteristics - The Skylab ATM CMG has a bandwidth 
approximately 30 radians per second when mounted on a rigid 
base. However, the bandwidth of the CMG is lower in this appli- 
cation since the CMG is mounted on shocks. Measurements made on 
the hybrid simulation show the CMG can be approximated as a second 
order lag with a damping ratio of 0.25 and a bandwidth of 20 rad/ 
sec over the frequency range of interest. More detailed analysis 
of the CMG is given in section 10. 
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Table 8-1. Closed Loop Gain at 0)=4ir for Various K^/K^ (Ideal GMG) 


b 

|T(jw)/T D (jo)) | 

20 log|T(jo))/T u (jo)) | 

0 

0.7071 

-3.01 db 

1 

0.7051 

-3.04 db 

2 

0. 7031 

-3 . 06 dib 

3 

0 . 7010 

-3.08 db 

4 

0.6990 

-3.11 db 

5 

0.6969 

-3.14 db 


Table 8-2. Factors of the Closed Loop Transfer for Various K^/K^ (Ideal CMG) 


b 

NUMERATOR 

0 

- 

s+4.3173 

1 

s+0. 24560 

s+4.0717 

2 

s+0. 52777 

s+3.7895 

3 

s+0. 87033 

s+3.4470 

4 

s+1.3464 

s+2.9709 

5 

s 2 +4.3!73s+5(£=. 

96538, o) n =2. 2361) 


s+1.2246 

s+1.6905 


DENOMINATOR 


s +8. 6346s+37.278(?=. 70711,0^=6.1056) 
s+0. 24515 s 2 +8 . 3895s+35 . 221 (C=. 70681 ,0^=5 .9347) 


s+0. 52271 s^+8.1119s+33. 

2 


=. 70565,w =5.7479) 
n 


s+0. 84363 s +7. 7910s+30.705(?=. 70300,0) =5.5412) 

n 


s +7. 4100S+28. 207 (?=. 69765,0) =5.3107) 
2 

s +6. 944ls+25. 539 (S=. 68704,0) =5.0536) 

n 





FREQUENCY - RAD 



PHASE - DEG 


Lgop Stability In vestigation - Stability analysis per- 
formed in this section is based on the assumption that the CMG can 
be approximated by a second-order lag with a natural frequency of 
13 rad/sec (handwidth-20 rad/sec) and a damping ratio of 0.25. In 
figure 8-1 the CMG is represented as C(s), the compensation network 
as G c (s), and the control law as G^Cs). 


The control analysis given in section 8.2 establishes the cri- 
teria for a system in which the CMG was assumed to have a transfer 
unction of 1. Since the shoekmounted CMG has the characteristics 
described above, a problem will exist which has not been considered 
xn the section 8.2 analysis. 


Writing the transfer function of the CMG as: 


C(s) = 




n 


s 2 -}-2C u s-Ro 2 
n n n 


, where : 


= 13 (rad/sec) 
C = 0.25 


( 8 - 22 ) 


the compensator as : 


G c (s)=l 

and the control law as: 

(8-23) 

s 

the open loop transfer function becomes: 


T GMG (S ) . (l+s/4.3) 
T (s) " 2 


[ 1+ 


0.5 ,s 
13 + ( 13 


J 


(8-24) 


The Bode diagram of equation 8-24 is given in figure 8-4. It is 
evident from the figure that the CMG has changed the stable system 
described in section 8.2.1 into an unstable system. Therefore 
compensation will be required to make the loop stable. Compensa- 
tion analysis is discussed in the following section. 
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Compensati on Analysis - Two candidates were considered 
' ’ The two were of the form: 


for the compensation network, G (s) . 


8-12 



GAIN (db) 



v s >=S! ^ 


(8-25) 


and 


G c ( S )= 


2C 1 S s 2 l 

1+ — — + -S-l 

03 2 

1 "lJ 



(8-26) 


TM c S h C ° nd ° rder com P ensator (equation 8-26) is used in this study. 
This type compensator is chosen so that the lead term of the com^ 

pensator can be chosen as the reciprocal of the second order lac 

L a ,^ eriStlC GMG ’ and the Gom P ensator lag can be 0.7 damp- 

ed with an of 100 radians per second. 

8 ! 6,1 Combined CMG/Compensator System - Since the CMG can be 
approximated by a second order lag over the critical frequency range 
and a compensator can be used as described in section 8.6, then the 
open loop transfer will again be described by equation 8-24. Now. 
however, S=0.707 and u> n =lO0 radians per second since the GMG lag 

and compensator lead cancel. Substituting these values into equa- 
tion a-24, the open loop transfer function is: 


G(s)=37.3 i 1+ . s /^-3) 
s 


1+.014s+1O _4 s 2 


) 


(8-27) 


gcneratJlY— Pl,>t ° f 8-27 and shous function 

60 degrees. 


* ~ ou \ j w» LUdL cms rune ci 

generates a gain margin (K ) of 24 db and a phase margin (0 ) of 

m 


A root locus plot is given in figure 8-6 of the rate plus posi- 
tion control law, the second order GMG and the compensator described 
by equation 8-26. The loop function for the root locus is: 


t cmg C s > 

V s ) 




r 2 q 
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(8-28) 
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where: 


C n 

= 0.25, 

(1) 

n 

= 13 

rad /sec 

? 1 

= 0.25, 

“l 

= 13 

rad /sec 

CM 

kj' 

= .707, 

U 2 

= 100 rad /sec 

R 

- V J 




a 

" V*r 





Rearranging equation 8-28 and substituting values suitable for 
constructing a root locus diagram: 


T fs) 

CMG KSJ „ ,s+4.3^, 
T(s) 2 K ' 2 


13 


s +6.5s+13 


2 2 
\ ^s +6.5s-hl3 x , 
K o > t 


13 


s 2 +140s+1Q0 2 


) C8-29) 


where : 


K, 


oj 2 r 


= 10,000 R 


Equation 8-29 shows the following pole and zero configuration 
exists. A double pole at the origin, poles at -3.25+J12.56, 
-7O.7+j70.7 and zeros at s=-4.3, s*=>-3.25+jl2.56. The one pole/ 
zero combination cancels, leaving the double pole at the origin 
and the zero at —4.3 essentially the same as described in figure 
8—2. The complex pole pair at s=— 70+j70 has a locus which crosses 
the jto axis and approaches the asymptotes at +60°. 

The locus near the origin will travel in a circular path. 
Therefore a value of gain (K^) will be chosen to put the closed 

loop roots on the 0=45° line to give 0.707 dominant root damping. 
This will give a set of closed loop roots at s=~4.6+j4.6 for a 
value of 1^=86,000. 


Using the assumption that the lead term of the compensator 
will nullify the CMG characteristics, equation 8-29 reduces to 


GCs) = 


i^Cs+4.3) 


2 2 2 
s (s + 2 c 2 w 2 s+t 0 o ) 


(8-30) 
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The closed loop transfer function of this 


compensated system 


is : 


t cmg ( s > 

T d (s) 


^(s+4.3) 


s 4 +25 2 oj 2 s 3 -H02S 2 +K ;l s4-4 . 3K X 


(8-31) 


Using previously chosen value of 1^, equation 8-31 becomes: 

(s+4.3) 


t cmg ( s) 

V S > 


4 3 * 2 " 

s +140s +10, 000s +86 , 000s+369 , 800 


(8-32) 


? f equatio11 8 ~ 32 factors into the four roots s=-4.6+14 6 
and s -66.1+J66.1. Therefore both sets of closed loop roots have 

^L“t„S1oi 0 - 707 - Equation 8 ' 33 che ^ ^ 


t cmg^ s ^ 

T d (s) 


K^s+4.3) 


2 ? ~ 

(s +9 . 2s+42 . 3) (s +l32s+8, 738) 


(8-33) 


A piotof equation 8-33 will give the system closed 
which is essentially the same as the solid lines on 
the frequency range of concern. 


loop response 
figure 8-3 over 


8,6>2 Combined System With Integral Gain - The 

< e PlUS law ualng 8 

.MG was discussed in section 8.2.2 Adding integral gain to the 
combined system discussed in section 8.6.1 will not affect the per- 
formance of the system appreciably. Equation 8-34 is the open loop 

law^with f U se tlG a a** 6 rate PlUS pos±t±on P lus integral control 
law with a second order approximation of the CMG and with the com- 
pensator described in section 8.6 


T 


CMG^ s) 

I £ (s) 


2 1 
s +4.3s+2 

r ^ i 
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2 2 -i 

s +25^8+0)^ 

r i 

i 
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a 

2 2 
s +2 Coo s+co 
n n 

2 

“l 

2 2 
Is +2 S s) w 7 s+(i) 7 . 


(8-34) 


Fi g ure 8-7 is a root locus diagram of equation 8-34. This diagram 
is essentially the same as figure 3-6, however, there are now three 

ZSS." 0 : L r c i6i Y nd “ Y 1 ' 1 ™ 1 — °» r«l 

°1« Jrinln ? ” lBl " “ d soes to °>« ™»1 near 

Scuf Y S - a " d Y Vel i n the left hal£ » la “ *» a Sear krc'Lr 

iocus. It is evident that the system would be unstable for a low 
or high value of ^ but if a value of ^= 86,000 is chosen^ the 

si!65?6i°65.6?° tS ^ f ° Und t0 ^ at S= "°* 525 ’ s— 4.05+j4.05, and 
The closed loop transfer function is: 


X CMG (S > 

t d (s) 


k^s+o^Hs+i.??) 


(s+0.525) (s 2 +8.1s+33) (s 2 +131s+8, 596) 


(8-35) 


07 S It ? ° f r °° tS in ec * uation 8-35 have a damping ratio of 
dotted ec luation 8-35 is essentially that given by the 

tted line in figure 8-3 over the 2 Hz bandwidth of interest. 

est , l-O-intini^rror Estimation - The pointing error can be 

estimated by assuming an ideal CMG. Considering the control law 

torn 5 P rT X ? , POS:Lt:ion (0) is 8 iven as a function of disturbing 
torque (T^) by equation 8-36. b 


0 (s) = 


Js 3 +K R s 2 +K pS +K I 


f T D (s)] 


(8-36) 


It has been established in section 5 that the major disturbances 
can be approximated by a step input with a magnitude of 15 n-m 
about the y and z axes. Applying the final value theorem 

lim 8 (t)=lim s 0 (s) 

t-x» S ->Q 

to equation 8-36 for a step input 


lim 0 (t)=lim{s 
t-w» 


s -0 Js^lC^+KjS+K 


db} 

s 


l±m 0 (t )=0 
t-x» 
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which says that the pointing error should be equal to zero for a 
step disturbance input when position integral is used in the con 
trol law. 

The preceding analysis is valid for linear systems, however 
there are nonlinearities present in the system which will cause 
pointing errors. These pointing errors will appear as limit 
cycles caused primarily by CMG friction. Results showing these 
effects are given in section 10 of this volume. 
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9. 


VEHICLE DYNAMIC MODEL 


pi 0r ff r to stud y the pointing performance realizable with the 

floated pallet under the control of the recommended CMC actuator 

necessary to formulate a model of the pallet/suspen- 
develonpH ^ yn ^ micS 7 t:o be used in conjunction with the CMC model 
invent-? t- a se ^ tlGn 7 ‘ s± nce only pointing performance is to be 
tudJ f’ he Vehlcle is essentially in an inertially held atti- 

tw'* a i? eref °T e 3 primary ^sumption used throughout this section is 

““dr” Vary Small> all °" lnB U ’ e S " aU 

Se parate rotational dynamics are developed for the pallet and 

forces'* d al ° ne W±th re a atiV£ translat ^ Disturbance to^^s and 
orces and gravity gradient torques act on both bodies, the CMG 

control torques act directly on the pallet, and forces and torques 

generated by the action of the suspension system act on both bodies. 

a The total slmula tion will include the digital CMG control law 

7 lnd P tL bffL” *•*! ° f tha GMG -«■ - given ^ section 
7, and the vehicle dynamic model of this section. The. pallet angular 

rates are picked off and processed through the vehicle control law 
^th^"^: ±n ° rd£r ^ “ torque commands for^se 

„ 9 ' 1 Resolu tion of CMG Control Torques - Although the control 

orque acting on the pallet is physically only an input to the dy- 

tJe Jndil'd i eDm p ined V6hicle * the of the torques from 

he individual actuators into a single torque vector is necessary, 

and conveniently discussed at this point. Given the output torque 
of the ith actuator as defined in equation 7-23 as T pI± , and using 
the CMG output to pallet^ transformations of equations 6-1 through 
6-4, the control torque ? cp acting on the pallet resolved into pal- 
let coordinates is: 


r CP = 2Z [ pVi J Vi 


i=l 


T = 
CP 


T P 'xl +T P«x2 +T P 'x3 +T P ' x4 

T P , yl +T P*z2' 1: p»y3 +T P' z 4 

T P * zl +T P 1 y2~ T P * z3~ T P * y4 


(9-1) 
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9,2 Derivation of t he Detailed Dynamics - The vehicle and sus- 
suspension system dynamics are derived in the following paragraphs, 
first with only the assumption of small angle restrictions in rota- 
tion and linear characteristics in the individual isolators of the 
suspension system. In section 9.3 the suspension dynamics are sim- 
plified based on the selection in volume III of an uncoupled four- 
point system. 

9 * 2,1 Suspension System Dynamics - When the pallet is in- 
ertially held under the control of the CMG system, the only connec- 
tion between the pallet and orbiter is through the suspension system. 
The suspension system exerts forces and torques on both the orbiter 
and pallet as a function of the relative rates and displacements of 
the two bodies. The dynamics of the suspension system are derived 
by considering each suspension point individually. It is assumed 
that each suspension point (isolator) exhibits only linear spring 
rate and damping characteristics along each vehicle axis with no 
coupling, i.e., a linear deflection or rate at any isolator loca- 
tion wili result in a force acting along the deflection or rate. 

In addition it is assumed that the individual isolators exhibit no 
rotational compliance or damping and finally that the undeflected 
isolator length is small compared to the distance from the pallet 
center-of-mass . 

In order to determine the forces applied to the orbiter and 
pallet by each isolator it is necessary to derive the deflection 

^1 of the ith isolator as a function of relative translation of the 

orbiter and pallet and the rotations of both bodies. Figure 9-1 
shows a representation of the isolator body geometry when the sus- 
pension system is relaxed and there is no translation of the orbiter 
or pallet with respect to the ideal (I) coordinate system. The 

vectors p pQ and p go locate the pallet and orbiter mass centers with 
respect to the origin of the I coordinate system. The vector R 
locates the pallet mass center with respect to the orbiter mass° 
center when the suspension system is relaxed. The vector I lo- 
cates the orbiter end of the ith isolator and a ± the pallet end. 

R o and ®i are fixed in fe he orbiter and rotate with it and a is 
fixed in the pallet. The vector P q coincides with R when the sus- 
pension is relaxed and is constant in I coordinates. 0 From the fig- 
ure, using assumption that the isolator relaxed length is negligible: 
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Figure 9-1. Geometry of One Suspension Point With No Deflection 


ITU ISOLATOR 
POINT 
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Now, consider the case when both the pallet and orbiter 
have been rotated with respect to 1 through the small angles 

0p ^ 0 Px ,9 Py ,0 p z ' ) and ^S = ^ 0 Sx’ 9 Sy’ 9 Sz^ ’ res P eciti vely , and both 
bodies transiated with respect to I coordinates as shown in fig- 
ure 9-1.^ The pallet and orbiter mass centers are located by the 

vectors and p^, respectively, with: 


P=P p -P s 


(9-3) 


The vector now represents the opening between the pallet and or- 

biter ends of the ith isolator, i.e., it is the deflection of the 
ith isolator. Expressing all vectors in their I space coordinates 

and noting that a., and R q have undergone rotations: 


P +I I T P 3 V [ I T S ] V l i 


(9-4) 


The vector £ is a measure of the translation of the pallet center 

ltS pos J- tion when neither orbiter rotation or relative 
orbiter/ pallet translation has occurred, thus: 


P=[ x T_JR+e 
Id o 


Substituting equation 9-5 into 9-4 and solving far t ! 

<V J t )+l i T pJ S i 

From equation 9-1, R -f =-a_. and: 

g i i 


(9-5) 


(9-6) 


Noting: 
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Figure 9-2. Geometry of One Suspension Point 
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defining: 

Equation 9-6 can be written as: 

-if -y 

o^e-hjjxa 


*i S <V*s>»i 


(9-7) 


(9-8) 


where is given in terms of its I space coordinates. Under the 
assumption of small angle rotations and small deflections £ the 
resolution of into orbiter coordinates is : 

i s = 1 s T i iVV¥A 

where 0 g x^ is essentially a second order terra and similarly for 

resolution into pallet coordinates. Thus, the isolator deflection 
can be expressed in both S and P coordinates by means of equation 
9-8, Under the small angle < assumption, ignoring second order pro- 
ducts, the deflection rate can be given as: 


• • • 

-f 

u^e+tjjxot^ 


(9-9) 
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The vectors and ^ are defined as directed from the orbiter 

to the pallet, resulting in forces on the orbiter along those vectors 
and in the same direction. Defining the spring constant and damping 
matrices describing the ith isolator characteristics as: 


tK i ] ' dla 8 fK *i-V K z i> 

(9-10) 

[c 1 ]^ia S {e jtl ,e y .,c zl ) 

(9-11) 


The force on the orbiter due to the ith suspension point is: 





(9-12) 


where the vector f i =( f xl ,f yi ,f 2i ) represents the friction of the 

suspension point, which for the purposes of this derivation is 
assumed to be of constant magnitude in each coordinate direction 
acting in the direction of the proper velocity component, i.e,: 


f k±= 5 ki >°' 

= 0 , 5 k± =0 • k=x,y,z 

= — F . , 5 , <0 j 

ki ki J 


(9-13) 


The force on the pallet due to the ith suspension point is simply 

the negative F„ . : 

So. 



C9-14) 


The torque exerted on the pallet by the isolator is due to the force 
-y 

acting on the lever arm a^: 


I pi = V F p i "- | V fF s 1 


(5-15) 


The torque applied to the shuttle is: 
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(9-16) 


T Si^i xF Si“ ^ R o + °i^Si 



R xF 
o 


SI 



The total forces and torques acting on each body due to n isolators 
as then: 


n 


SS r Si 

(9-17) 

i=l 


F =-F 
SP SS 

(9-18) 

n 


SP 1 Pi 

(9-19) 

i=l 


ss = \ xF ss“ t sp 

(9-20) 


where the subscript "SS" indicates action of the suspension on the 
orbiter and "SP" on the pallet. 

9 ' 2 ' 2 Qtbiter R otational Dynamics - Under the assumption of 
small angle rotations, the rotational equations of motion of the 
orbiter are simply: 




(9-21) 


where the Euler term involving the cross product of rates has been 
neglected, and the external torques are the suspension torque T 
(equation 9-20), gravity gradient torque acting on the orbiter 
^GS an< ^ disturbance torque acting on the orbiter T^ g . The matrix 
[J g ] represents the inertia tensor of the orbiter. The orbiter 
rates and displacements t g and $ g are obtained by integration of 
the acceleration components. 
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“pHVA (9-22 

where [J p ] represents the Inertia of the pallet and the external 

and^eontrol gradient, ^turbanee 

Pallet rates a„ d P dlspIa’L^ m 

or equation 9-22. rained from the time integrals 

«- Pailat JJ. arbiter and 

complete solution of the translational t* 1 coordlnate frame, a 

"he relative ' transit ^ RS 1 T 

bh. arbiter is retired to £ 

particular the vector e she™ In figure 9-2. In general, the vector 
e expressed in I spaC e components can be related to p and R , re- 
calling R g rotates with the orbiter by: °* 


P^E+fj-TgjR -e+R +0 XR 
r o o o S o 


from which; 


-y -y 
£=p-R -I xR 

o S o 


and, recalling R q is a constant magnitude 
the small angle assumption: 


vector, and again using 


e=p-9xR 

s o 


(9-24) 


It should be noted here that the vector ? c which appears on figure 
8-2 which is given by V«V\ is «“ ^ive translation used 

the two derivations. compatibility between 
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- Aea±n ’ r ' ferrl "8 <=» «S“» 9-2 it tan ba seen that tha vector 
P can be expressed in terms of J and J which locate the pallet 
and orbiter mass centers as: 


P=P p -P s 


and thus , 

~r “>■ •> 

p=p p -p s 

Considering each body separately, the 
iation are: 


(9-25) 

equations of motion in trans- 


Vp*V ? DP 

Hj -► 

m„ p =F +F 
S H S SS DS 


(9-26) 

(9-27) 


where m, and m s are the masses of the pallet and orbiter, F and 

P s p the forces of the suspension and * „ p and ^ disturbancfforces. 

olving equations 9-26 and 9-27 for the acceleration, substituting 
into equation 9-25 recalling F =-F 


‘SP ~*SS - 

g- F pp _ f ds _ * 

“p S Vs Fss 


(9-28) 


The time integral of equation 9 7 28 gives g for use in equation 9-24 
since the initial condition on p=0. However, to determine p: 

t • 

P" J P(T)dx+p 
0 G 

and, since tha I space representation of S is the same as the s 


space representation of R ; 

o 


J P (t) dT=>p-p Q £e^=p_p 


(9-29) 
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equation 9-23 to deter- 


Thus, the integral of p Is substituted into 
mine e. 


°rqua; 2 ;ctlnr^ r ^hf - AH external forces and 

-orques acting on the bodies are put into the equations of motion. The 

pallet rates 0 p are picked off and processed through the vehicle con- 
trol law, GMG control law and CMC dynamics to obtain the control torque 
T cp . The suspension forces and torques are fed back into the dynamics 

sznzzs: In vector 


Pallet Rotation: 


(9-30) 

Orbiter Rotation: 

V[ J s rl ( J SS +f GS +:f D s) 

(9-31) 

Translation: 

“p^s t 

“p m s ®p® s SS 

(9-32) 

Integrating: 

V V T %-/ t p dT 

(9-33) 


e s -i ^s dT t s -l ® s dT 

(9-34) 


p=/ PdT e=p-R =/ PdT 

Li 0 

(9-35) 

The suspension dynamics are then determined by: 



HA H 4 s 

(9-36a) 


• • « 

^ ■** ->- -v -v 

e-e c 4 s xR o 

(9-36b) 

For i=l, n where n 

is the number of suspension points : 
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^ 1 =£-+^xa 1 

? si =fK i ] V tc i j V ? i 


T Pi~ _a i xF Si 


(9-37) 

(9-38) 

(9-39) 


and torques to be fed back to *1“““°™ 


n 


F = > F 
SS 2-_. Si 


z 

i=l 

n 

sp Z_, x pi 
1=1 


(9-40) 


(9-41) 


VW*sp 


(9-42) 


, n , 9 ' 3 „ Z ^plified Veh icle D ynamics - The previous derivation if 

woul’d me T 1 l d f ° r s±TOulation would result in a dynamic simulation which 
would allow parametric variation of the number of isolators and t-hp-fr 
characteristics and locations, as well as canLrfof "a“ of 

oth the orbiter and pallet. Within the scope of this study however 
such generalizations are not needed and significant savings can be ’ 

and'uJinB thrfo U£atl i n t reI,Ulred by makln 8 additional assumptions 

thin: P i “t: uspBnsi °” * 

q o ^e four-point configuration is shown schematically in figure 
-3 which illustrates the near symmetrical placement o/the iso- 
lators about the pallet center-of-raass and in a plane containing 
the center-of-mass. Under these conditions, the locations of Se 

isolators a ± , ±=1,4 can be expressed in terms of two lengths, i.e., 

the longitudinal and lateral distances from the mass center ’ As-" 
uming perfect symmetry which implies uncoupled translation and 

rotation for the pallet, the isolator locating vectors « can be 
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"p 



Figure 9-3. Suspension Point Identification - Four Point System 
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given in terns of the location of the number 1 isolator. Notion 
he isolators are numbered in terms of the quadrant of the /y 8 

Plane in which they lie, the S ± become under this assumption] P 




«T<W 0) 

(9-43a) 

V < "VV 0) 

(9-43b) 

V ( “V"V 0) 

(9-43c) 

V (< V" a y' 0) 

(9-43d) 

4 


i=l 

(9-44) 


[C ± ] [C] and F k .=F k , the expression for the forces and torques ap- 

as sumptions 6 '^ P ? nSiGn Can be 8 reatl y simplified. With these 
assumptions, the force equation 9-40 can be written: 


F ss = {TK-JCe^x^J+lcJCe+Ixa^+f } 

i a l 

4 4 4 

“4 [K]e+[K]$x a^+4 [G] e+[g] ^ a ± + ^ J 


or, 


F ss =4[K]e+4[Cjt+ ^ ? 


(9-45) 


It should |e noted that although the quantities 6 are no longer 
required, 5 ± must be calculated to determine the frictional forces 
f i, however, expanding the various $ : 
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-a 

y z 

V 

a 4 
x z 

• • 

„ a (p -a 4 
[ y x x Y y 


^ 2 =e+ 


-a 4 
y z 

-GX 4 
X z 

• a 

a <j> +a 4 
x y y x 


1— E+ 


a 4 
y 2 

-a 4 

X z 
• • 

a <J> -a 4 

x y y x 


, = £+ 


a 4 
y z 

• 

a (b 
x z 
• • 

-a <p -a d> 
x y y x 




By making the following definitions of individual components: 


i =e -a i 
xl x y z 

• • • 

o =£ +a 4 
x3 x y r z 

6 =e +a 4 

yl y x T z 

• • • 

q =£ -a 4 
y2 y x z 

• • • a 

o ,=£ +(ct (j) -a 4 ) 
zl z y x x Y y 


(9-46a) 

(9-46b) 

(9-46c) 

(9-46d) 

(9-46e) 
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(9-46f ) 


6 z2 =£ z + ^v^ x V 

*z3 E z ^°y^x a x^y) (9-46g) 

6 z4 _e z~ (a y^x +a x <l> y ) (9-46h) 

The twelve individual components of the <$. can thus be given in terms 

ffrence a of V lnv 8 li G ^ lGU ^ tlons which in turn involve sums and dif- 
many identical quantities. Denoting f as the friction 

resulting from the rate <5 , a=x,v z* i=l 2 A ►u,, 

ai ,y,z > 1 -*->^>->>4, the total friction 

force becomes : 


t t * E t- 

f i=i 1 


2(f xl +f x3> 

2 <W 

f zl +f z2 +f z3 +f z4 


(9-47) 


Consider next the torque equation 9—41: 


C SP ^ a i x ^ M (e+4>xa)+[C] (E+|xa )+f } 
, . "*■ 1. 1 


=[K]ex 


^ V [K ^^±> + fC3ex a ± - ^ a ± x[c ](fcS )- Y Z.xf 

1 11 1 1 


The terms involving only sums of the a ± clearly disappear, then ex- 
panding : 


2^ a^xlK] (|xa i )= Y 


. K (a <j> -a ,d> ) 

4 z vi x xi 


z yi r x xi~yiV 

2 

K (a tj> -a ,a d> ) 

z xi Y y xi yi^z' 

(K a 2 +K a 2 )<b 

x yi y xi'^z j 


4K a 2 d> 

z y r x 

4K a Z (j) 
z x y 

4(K cx 2 +k: a 2 )(b 
I x y y x ' v z 
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“io" ?o?"e= lel £or «« <»«**■* te™. Defining . fric . 


m Asr^-»- -+• 


Then, defining: 


- a y (f Z l +f z2- f z 3- f z 4) 

a x (f Z l- f 2 2- f z3 +f z4> 

L 2a y ( f xl" f x3 } “ 20t x (f yl -f y 2 > 


[K c J=4[KJ 

£D j^4fGj 


f K (j) J=diag{4a 2 K ,4a 2 K 4(cx 2 K -hA )} 
A i . I . J X X 7 


..2„ . 2 
x i 

(B^diag{4a Z B 4aJ 4 (a 2 D -kx 2 D )} 

T y * x z yxxy 

The suspension torque and force equations become: 

VV^VH 


T sp -- [K<^ J ] |-T f 


and as before. 


(9-48) 


'9-49a) 
(9-4 9b) 
(9-49e) 
C9-49d) 

(9-50) 

(9-51) 


T SS" T SP +R o xF S 


(9-52) 


The equations of motion cannot be further simolif^H -in 

however, since the cross products of inertia are ^ 

Sfs^plSefd^c 6 IT" 8 3 dia ^^’^a G tenLr 
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Orbiter / Suspension/Pallet Dynamics 
Figure 9-4. Vehicle Dynamics Block Diagram (Sheet 1 of 2) 




































10. COMPUTER STUDY 


This section, gives the computer results obtained from the 
hybrid simulation of the orbiter vehicle and the suspended pal- 
let controlled by the shockmounted ATM CMC. The results described 
in this section were obtained from the pointing accuracy model 
defined in the preceding sections with the exceptions that CMC 
gimbal angle change, translational dynamics and suspension torque 
on the shuttle were neglected since they had negligible effect 
on the overall system response. 

The main purpose of this study was to determine pallet 
pointing capability as a function of CMC nonlinearities, crew 
motion disturbance, GG effects, and system bandwidth. Studies 
determined that the only significant CMG nonlinearity was fric- 
tion, discussed in section 7,3, and that GG disturbances would 
not affect pallet pointing appreciably. Therefore pointing 
studies were parameterized as a function of friction level, 
crew motion disturbance, system bandwidth and CMG gimbal angle. 
Limit cycle levels were obtained with no torque input and then 
the torque input level required to eliminate limit cycles was 
determined. These results are presented in this section with a 
summary of results given in table 10-4. 

10.1 Hybrid Simulation Details - Figure 10-1 is a block 
diagram of the complete simulation showing the interrelation- 
ships of the vehicle control law (equation 8-1) , the compensa- 
tor (equation 8-26), the CMG control law (equations 6-13, 6-14, 
6-22, 6-24 and 6-25), torque resolution (equation 9-1), the 
vehicle dynamics and the CMG models, both of the latter are dis- 
cussed below. The model was implemented on a hybrid computer 
facility comprised of a Scientific Data Systems 86 digital com- 
puter interfaced with two Applied Dynamics AD-4 analog consoles 
and a Beckman Ease analog computer. The simulation was run in 
real time with the fastest possible sample period which was 

7 msec. No quantization was present in the runs other than 
that of the computer A/D and D/A interfaces, approximately 0,006 
percent of maximum value. 

10.1.1 CMG Models — A detailed block diagram of the CMG 
model is given in figure 10-2. The only simplification of the 
model derived in section 7 is the deletion of the various limi- 
ters present in the rate loops. This was done since the region 
of operation of the CMGs in this study involves low signal levels 
throughout, thus obviating the need for the limiters. The CMG 
parameter values used in the frequency response runs and in all 
pointing performance runs with the exception of those involving 
variations of the CMG friction are listed in table 10-1, 
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tlon were^T 1 ^^ m ° del * inc °rporated into the complete simula- 
tion were altered somewhat to minimize the analog hardware re- 
quirement. This simplification consisted of replacement of the 

IttLTjt r 5 ° f 6l,tbal a " 8leS “ ith «“»« gains which could 

ficaMn re P resent any arbitrary gimbal state. This simpli- 

eimhln r S % P °f iBle due to short time periods and slow 
gimbal rates involved which validates the use of linear approxi- 
mations to these functions. approx^ 

„„ 10,1 ; 2 Vehicle Dyn amics Model - Figure 10-3 is a block 
diagram of the simplified vehicle dynamics used in the Simula- 

SL'd y ' a The * y ™'° 1CS tln * lly ln s“dy arc 

those derived in section 9.3 with the translational dynamic! and 

suspension torque acting on the crbiter being neglected due to 
the very low magnitudes of these items compared Jith the other 
parameters existing in the simulation. 

in !°~ 2 1±S 5 S thE vehicle dynamics parameter values used 

xn all pointing performance runs. 

- 10,2 . - MG Frequency Response - The CMG model frequency response 

fol^r X T r f d ° uter simbal angles and 0.01 rad/sec gilbal rate 

T ^rthe ra ™? tU ^ g±Ven in figures 10 “^ and 10-5. The quantity 
c the CMG output torque in the direction of the torque command 

tnL P o7^ diCUlar f thE CMG m ° ,nentum vector divided by the ampli- 
tude of the command torque and is theoretically one. The quantity 

x is the CMG output torque in the direction perpendicular to both 

the command torque and the CMG momentum vector divided by the am- 
plitude of the command torque and is theoretically zero. The third 
rt? mP rMp nt ° CMG output Corc iue, the component in the direction of 
monitored!^ 11 " 1 W “ J ' udsBd t0 be 3raa11 and therefore not 

tho CGmp£nsate for the CMG dynamics of figures 10-4 and 10- 5 
the following compensator parameters (equation 8~2g) we-e -mp loved 
^n all pointing performance runs: 10^=12.5 rad/sec, £ =0.125^ 7 

rad/sec and C 2 = O.707. 1 

widthf;' 17 QUtllna “ Pallet £ ° ntrGl loGp ba'nd- 
xacns ot i/2, 1 , 2 and 4 Hz were considered in this stndv tv, q 

T tI01 la ! V° SU1 ™- P-iti"7J,”egrS y sal L 

c o r r e s p o nd ine to these bandwi rlfr-fac Mam • i r 


, . ” - — > ana position integral erains 

. ^responding to these bandwidths were determined from the^naly- 
sis of section 8 and are given in table 10-3, y 
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bandwidths an^at^cm eimbal° f made f ° r th ® varlous system 

T ha M aC UIG Smjbal angle combinations of 0° and 

The_th.ee types of pointing performance runs which were mad f * 

1" W ° tQrque dis turbance acting on the pallet or arbiter. 

2. Step torque disturbance acting on the pallet. 

tn Jr? ^ iS 1 iHdlviduall y at the minimum value necessary 
to eliminate the undisturbed limit cycle on that axis ! 

pirt^! aXeS simultaneous ly at the values determined in 

3 ; GreW -^torque disturbance acting on the orbiter as 
given by where R=(-15,0,0) T meter, and 

Case 1 v^v°) T n 

Case II F ds =( 0 ,0,F d ) T n 

Gasem ^=( 0; |^ T n 

where: ^ ramps from zero to 100 newtons during the first 

0.8 seconds, is zero between 0.8 seconds and 24 
seconds ramps from -100 newtons to zero during 
the las t 0 . 8 seconds * ■ 

The above runs were made for: 

a) The 1/2, 1, 2, and 4 Hz systems with all gimbal 
angles equal to 0°. gimoai 

eq U al h to 1 M°? , ’ d 1 ^ SySt ™ S " ith 311 8lmbal a "Sles 

equal h to 2 0””f„K H % syst f"f » lth the inner gimbal angles 
qual to 0 and outer gimbal angles equal to 45°. 

systems T**/™ fhe 1, 2 and 4 He 

adjusted to 1/2 the nominal friction le“l Sd to 2 Slf^e"” 
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nominal friction level. Time histories of various system vari- 
ables were plotted on strip chart recorders. The variables 


chosen for display were 

the following: 

y j 

= 1,2, 3,4 

CMG inner gimbal rates 


= 1,2, 3,4 

CMG outer gimbal rates 

^Si’ 1 

= x,y,z 

orbiter angular rates 

0 si ji 

= x,y,z 

orbiter attitude 

T DPi ,;L 

= x,y,z 

torque disturbance on pallet 

T Ci ,:L 

= x,y,z 

torque command from vehicle control law 

’Si ’ 1 

= x,y,z 

torque command from compensator 

T CPi* 1 

= x,y,z 

net CMG torque applied to pallet 

^i ’ 1 

= x,y,z 

pallet angular rates 

Vi ’ 1 

= x,y,z 

pallet attitude 


10.4 Pointing Performance Study Results - Reproductions of 
selected computer runs are shown in figures 10-6 through 10-46. 

In. reading the runs, time can be read from the marks at the top 
and bottom of each run where each mark corresponds to 1 second. 

10.4.1 CMC Friction Study - The primary nonlinearity of 
the CMG is friction. When no vehicle disturbances are present, 
the CMG friction causes limit cycles as evidenced by observing 
Q pi on figures where no disturbance torque is applied (i.e. , 

figure 10—6 for example) . As can be seen by comparison (figure 
10-6 vs 10-16), the limit cycle amplitude decreases as the system 
bandwidth increases. Figure 10—47 summarizes how pallet attitude 
(0 pi ) varies with bandwidth. 

Studies were also made to determine the effect of friction 
level on limit cycle amplitude, (See figures 10-41 through 10-46.) 
For a 2 Iiz system with gimbal angles set equal to zero, the pallet 
attitude error due to nominal CMG friction was approximately 0.1 
arc-second about the pallet Y and Z axes, and approximately 1.1 
arc-second about the pallet X-axis. As the CMG friction level 
increases or decreases, the limit cycle amplitude varies propor- 
tionately. The effect of friction level for various system band- 
widths is summarized in figure 10-48. 
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torques ' fT 2 ^^bance ~ When snail vehicle disturbance 

q S ^DPi 3 are a PP lled > approxiamtely 1 n-m, steady CMG gimbal 

rates occur and the system stops limit cycling (i.e f inure in 
Disturbance torques (T^) required to stop limit ciciiJS”^^ 

system bandwidths are summarized in figure 10-49, 

. 1G * 4 * 3 — ew Mqtion Disturbances - Crew motion is the lamest 
disturbance expected which will affect pallet pointing performance 
the 2 Hz system with CMG gimbal angles at zero, the worst case 
crew motion disturbance error is 0.25 arc-second about the pallet 

pointing 366 As can be seen from the computer runs 

err ° r lnCr °duced as a result of crew motion disturbance ’ 

wt^rLIlgur'S'o-M lt,Cr ° aS “- These are sum- 

3-0.5 S ummary - Data presented in this section shows that the 
system will have a limit cycle attitude error due to the CMG fric- 
ion when no external torques act on the vehicle. It also shows 
that the limit cycle amplitude to be a function of friction level 
and system bandwidth. This limit cycling can be stopped by apply- 

required UeS In°add^tT ehiGle Wh±Ch CaUSe CMG girabal raCes to be 
emS rhich i r\ CreW m ° tl0n causes attitude 

error which is reduced by increasing system bandwidth. 

seem that^h^h* h** Sa±d ±n thG above Paragraph, it would 

seem that the higher the system bandwidth the better the system 

ESrsr h r e T* the higher tue nt h ii iiit 

at 45 delree 7 ? ^ In fact « the 2 Hz system was unstable 

aS 5 ae lnn f r and QUter eimbal angles and the 4 Hz system 
was unstable ^ at both 45 degree inner and outer gimbal angles and 
zero inner gimbal angles and 45 degree outer gimbal angles. This 
instability was not predicted by an idealized CMG system a^lysL 
and thus must be attributed to CMG dynamics. However, it is anti- 

SfSectfo? thfcMf f f 6 CMG l0 ° PS> taklnS int ° aacount 

ties shock-mount, will eliminate these instabili- 


su mmary, the 2 Hz system could easily meet the 1 arc- 
second pointing requirement considering the fact that an orbital 
vehicle almost always has disturbance torques acting upon It 
Redesign of the CMG rate loop is reco^JL t “taS S amount 
. S effeGts bb e CMG shoekmounts and eliminate system instabili- 
ties presently observed for the ATM CMG rate loop design. In^ddi- 
tion a redesign of the rate loops, to better compensate for the 
fects of friction at zero gimbal rate commands, will allow lower 
vehicle loop bandwidths (i.e., approximately 1 Hz) while still 

srss s:tL~? nd pointin8 stab±iity 
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Table 10-1. Floated Pallet Pointing Performance 
Study CMG Parameters 


Parameter 

Symbol 

Nominal 
Magni tude 

Units 

Inner Axis Compensator Gain 

K 1 

10.54 

vde/vdc 

Outer Axis Compensator Gain 

K 3 

12.80 

vde/vdc 

Inner Gimbal Amplifier Gain 

K al 

10.0 

amp/vde 

Outer Gimbal Amplifier Gain 

K a3 

10.0 

amp/vdc 

Inner Gimbal Motor Gain 

K ml 

1.559 

N-m/ amp 

Outer Gimbal Motor Gain 

K m3 

1.559 

N-m/amp 

Inner Gimbal Tachometer Gain 

K R1 

0.995 

vdc/rad/sec 

Outer Gimbal Tachometer Gain 

^R3 

0.995 

vde/rad/sec 

Compensator Lead Corner Frequency 

“n 

20.0 

rad/sec 

Compensator Lag Corner Frequency 

“d 

0.2 

rad/sec 

Motor Deadzone 


0.003 

ampere 

Gimbal Motor Static Friction 

s 

in 

0.1 

N-m 

Gimbal Motor Running Friction 

R 

0.08 

N-m 

Zero Gimbal Rate Band 

E g 
E m 

0.00005 

rad/sec 

Zero Motor Rate Band 

0.002828 

rad/ sec 

Gear Ratio 

N 

g 

^ml 

56.55 


Inner Gimbal Motor Inertia 

0.0068 

, 2 
kg-m 

Outer Gimbal Motor Inertia 

J m3 

0.0068 

. 2 
kg-m 

Inner Gimbal Gear Train Stiffness 

K gl 

203,370 

N-m/rad 

Outer Gimbal Gear Train Stiffness 

K S3 

L 

203,370 

N-m/rad 

Moment of Inertia of Wheel About 
Axis Perpendicular to Spin Axis 

1.803 

kg-m 

Gimbal Static Friction 

S 

g 

R 

R 

0.5 

N-m 

Gimbal Running Friction 

0.4 

N-m 

Bearing Compliance 

K 1 
yi 

%3 

542,320 

N-m/rad 

Bearing Compliance 

542,320 

N-m/rad 

Compliance Damping 

v 

20 

(N-m-see) /rad 

Compliance Damping 

6 Y3 

20 

(N-m-sec) /rad 

Angular Momentum of Wheel 

H 

r 

3115.0 

N-m-sec 

Inner Gimbal Inertia 

J all 

0.746 

2 

kg-m 

Inner Gimbal Inertia 

J a22 

0.868 

2 

kg-m 

Inner Gimbal Inertia 

J a33 

1.003 

. 2 
kg-m 
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Table 10-1. Floated Pallet Pointing Performance 
Study CMG Parameters (Concluded) 


i 

i 

1 


Parameter 


Outer Gimbal Inertia 

Outer Gimbal Inertia 

Outer Gimbal Inertia 

Equivalent Outer Gimbal 
Inertia 

CMG Base Inertia 
CMG Base Inertia 
CMG Base Inertia 


Shockmount Spring Const. 
Shockmount Spring Const. 
Shockmount Spring Const. 
Shockmount Damping 
Shockmount Damping 
Shockmount Damping 


Symbol 

Nominal 

Magnitude 

Units 

J ell • 

0.76 

, 2 
kg-m 

J c22 

4.3 

. 2 
kg-m 

J c33 

3.9 

. 2 
kg-m 

J c33 

4.834 

kg-m 

J bll 

5.35 

kg-m 

J b22 

7,9 

. 2 
kg-m 

J b33 

3.63 

t 2 

kg-m 

K Sx 

170,000 

N-m/rad 

K Sy 

208,000 

N-m/rad 

K Sz 

58,000 

N-m/rad 

D Sx 

402 

(N— m-sec) / rad 

D Sy 

490 

(N— m-sec) /rad 

D Sz 

137 

(N-m-sec) /rad 
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Table 10— 2. Floated Pallet Pointing Performance 
Study Vehicle Parameters 


Parameter 

Reciprocal of Pallet Inertia 
Reciprocal of Pallet Inertia 
Reciprocal of Pallet Inertia 

Reciprocal of Shuttle Inertia 
Reciprocal of Shuttle Inertia 
Reciprocal of Shuttle Inertia 

Suspension Spring Const. 
Suspension Spring Const. 
Suspension Spring Const. 

Suspension Damping 
Suspension Damping 
Suspension Damping 


Symbol 

Nominal 

Magnitude 

A Pxx 

4.696 x 10" 5 

A Pyy 

7.204 x 10~ 6 

^zz 

7.384 x 10“ 6 

A Sxx 

1.015 x 10~ 6 

A Syy 

1.385 x 10" 7 

A Szz 

1.354 x 10~ 7 

K 

8,420 

K, 

d>y 

21,555 

K. 

9 2 

53,432 

c 4>x 

2,679 

c 

4>y 

6,860 

c, 

<pz 

17,007 


Units 

1/kg-m 2 

l/kg-m 2 

l/kg--m 2 

1/kg-m 2 

1/kg-m 2 

1/kg-m 2 

N-m 

N-m 

N-m 

N-m-sec 

N-m-see 

N-m-see 





Table 10-3. Floated Pallet Pointing Performance Study 
Vehicle Control Lav: Gains 
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Table 10-4. Floated Pallet Hybrid Simulation Study Results 




N N 
fd « 


T3 

s § 


w 

H C3 

£9 

3*(S 


95 


(PEAKVALUE) 

SEC 


Case III 
(Y axis) 
(Z axis) 


SYSTEM BANDWIDTH 0.5 Hz 
"1 

1.0 Hz 

2.0 Hz 

4.0 Hz 

2 TIMES 
NOMINAL 
FRICTION 

X axis - 

Y axis 

Z axis - 

mm 

i 

1 

mm 

mM 

ib^b 


X axis 4 . 0 

1.5 

■ 

■ 

NOMINAL 



i 


FRICTION 

Y axis 0 . 4 

0.15 


0.05 


Z axis 0.4 

0.15 

0.1 

0.05 



0.5 TIMES 

X axis 

- . 

0.7 

iwm 

0.4 


NOMINAL 

FRICTION 

Y axis 

- 

0.1 

0.05 

0.03 



Z axis 

- 

0.1 

0.05 

HRS 

TORQUE DISTURBANCE 

X axis 

0.5 

0.6 

1 fi 

2.2 

REQUIRED TO 

Y axis 




ELIMINATE LIMIT 

1.2 

1.1 

0.8 

1.4 

CYCLE 

(N-m) 

Z axis 

1.2 

mm 

0.6 

1.4 



Case I 

m 




CREW MOTION' 

(Z axis; 

0.25 

0.05 

DISTURBANCE 

Case II 

B 




(PEAX_VALUE) 

(Y axis) 



nil 

SEC 














(X axis) 

■SB 


0.05 

nil 

“ ■ . . . — — 


(Z axis) 

IK&fl 


0.1 

0.05 

UNLOADED LIMIT 

X axis 

im 


0.5 

A 

CYCLE 






(PEAK VALUE) 

Y axis 

0.5 

0.25 

0.2 

* 

SEC 








Z axis 


0.2 

■n 

* 

TORQUE DISTURBANCE 

X axis 

5.0 

mm 



REQUIRED TO 


IBB 

"k 

ELIMINATE LIMIT 

Y axis 

1.2 


1.2 

k 

CYCLE 

(N-m) 








Z axis 

1.2 


1.0 

* 



Case I 
(Z axis) 

5.0 

Dii 

0.2 

* 

CREW MOTION 
DISTURBANCE 

Case II 
(Y axis) 

2.0 

0.5 

nil 

* 


nil 

0.15 


Unstable system, not run. 
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Figure 10-1 - Floated Pallet Pointing Performance Study, 
Hybrid Computer Model 
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Figure 10-3 - Floated Pallet Pointing 
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Figure 10-12. Step Torque Disturbance, 
45 deg Gimbal Angles 
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Figure 10-14, Crew Motion Torque Disturbance Case II 
1/2 Hz System, 45 deg Gimbal Angles * 
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Figure 10-15. Crew Motion Torque Disturbance Case III 
1/2 Hz System, 45 deg Gimbal Angles 
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Figure 10-26 . No Torque Disturbance, 2 Hz System, Zero Glmbal 


Angles 
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Figure 10-34. Crew Motion Torque Disturbance Case II, 2 Hz System 

Zero Inner Gimbal Angles and 45 deg Outer Gimbal Aneles 
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Figure 10 50. Crew Motion Disturbance Peak Value Versus 


System Bandwidth (Zero Gimbal Angles) 



11. CONCLUSIONS 

Ihe 1 ollowing are the major conclusions of this study: 

a. Two shuttle pointing orientations (X-IOP and X-POP) were 
studied to determine actuator torque and momentum storage 
requirements. The X-IOP attitude has the advantage of 
allowing instrument pointing anywhere in the celestial 
sphere without the need for gimballed payloads, however, 
momentum accumulation due to gravity gradient torques 

is larger than when operating in the X-POP attitude. 

It was derided that the momentum exchange system would 
oe sized for X-IOP in order to allow all sky viewing 
without the requirement of additional gimbals when con- 
sidering single payload operation. 

b. A CMC cluster having a momentum storage capability of 
12,200 n-m— sec and a minimum torque capability of 200 
n-m can meet the requirements of the orbiter /pallet 
vehicle for gravity gradient stabilization during obser- 
vation in an X-IOP orientation, and for momentum desatura- 
tion maneuvers. 

c. Of the three candidate double gimbal CMC actuators 
(Bendix MA-2300, Bendix MA-2000, Sperry 4500) the modi- 
fied Bendix MA-2300, a version of the Skylab ATM CMG was 
chosen primarily on tho basis of being space qualified 
and its off-the-shelf availability. This choice dictated 
a CMC cluster of four Skylab ATM CMGs of 3,li5 n-m-sec 
momentum capability each in order to meet vehicle momentum 
requirements. 

d. A 20 Hz isolator was placed between each CMG and the 
pallet. High frequency disturbances generated by the 
CMG, which could not be eliminated by a second level of 
Experiment Control (e.g., Image Motion Compensation) due 
to its high frequency nature, are prevented from reaching 
the pallet. 

An unwanted effect of the isolator is to lower the 
frequency response of the actuator (CMG/isolator combina- 
tion) to approximately 2 Hz. Design of a compensator was 
required to insure stability when operating at a system 
bandwidth of 2 Hz. 
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e. 


f. 


Computer studies were parameterized as a function of CMr 
friction level, CMC gimbal angle, syste „ JF 

crew motion disturbance. Major conclusions obtained from 
the computer runs were: 1 ) with no torque present there 

exists a pallet attitude limit cycle Hi JL a > 

mum for 2 Hz system) which decreases in'amplUude' vUh“„ 

varleH 36 “ SyStem band width, 2) limit cycle amplitude 
varied propo.tronateiy with friction level, 3) a torque 

tifi 1 !" Ce ( f p P roximate ly 1 n-m for a 2 Hz system) causes 
the limit cycle to disappear, the torque level being a 

function of system bandwidth, and 4) pointing error due 
*° " 1 ° tl0n dis turbance (0.25 arc-second for a 2 Hz 

y em) decreases as system bandwidth increases. 

theMK r f T rate l00ps ’ t0 better compensate for 
effects of the CMG shockmounts and friction at zero 

gimbal rate command o, will allow the use of lower vehicle 
ine P che n ? WldthS (i ' e ’ ap P roxilna tely 1 Hz) while still meet- 
for all iys a t:m S :ta°tesr° lntin8 StjbUity perf «e desired 
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12. NOTES 


12 . 
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2 . 


3. 


4. 


5. 


12.2 

A/D 

APCS 

ASM 

ATM 

CM 

C&D 

CMG 

D/A 

DCMG 

GG 

kg 

kg-m' 

km 

n-m 

n*mi 

OMS 

RAM 

RCS 
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Abbreviations 


Analog to Digital 

Attitude Pointing Control System 

Astronomy Sortie Mission 

Apollo Telescope Mount 

Center-of-Mass 

Control and Display 

Control Moment Gyro 

Digital to Analog 

Double Gimbal Control Moment Gyro 

Gravity Gradient 

Kilogram 

2 

Kilogram- (Meters) 

Kilometer 
Newton-Meter 
Nautical Mile 

Orbital Maneuvering Subsystem 
Research and Applications Module 
Reaction Control System 
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r Pm Revolutions Per Minute 

SCMG Single Gimbal Control Moment Gyro 

SRB Solid Rocket Booster 

wrt With Respect To 

12.3 Symbols 


a 


i 


A~space 

B-space 

C-space 


[CJ 

C(s) 

D si 


-y 

F„ 



F(s) 



G 1 (S) 

G c (s) 


Components of unit vector along the local 
vertical (i=x,y,z) 

CMG inner gimbal coordinate frame 
CMG base coordinate frame 
CMG outer gimbal coordinate frame 
Pallet isolator spring damping matrix 

CMG/shockmount transfer function 
Shockmount spring damping 

Total friction force on orbiter from suspen- 
sion system 

Friction of the ith pallet suspension point 

Force on the pallet due to the ith suspension 
point 

Lapace of any time function f(t) 

Force on the orbiter due to ith suspension 
point 

Sum of the suspension forces acting on the 
pallet 

Sum of the suspension forces acting on the 
orbiter 

Vehicle control law 
Compensation network 


H 


Aj 


-y 

H 


—y 

H. 


H 


gg 


-y 

H 


ggm 


Angular momentum of the inner gimbal written 
in A-space 

Wheel momentum of the jth CMG (i=l,2,3,4) 

Angular momentum of the outer gimbal written 
in C-space 

Accumulated momentum for one orbit 

Momentum accumulation due to gravity gradient 

Upper bound on momentum due to gravity gradient 
torque 

Peak momentum due to gravity gradient torque 
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-V 

H 

r 


Angular momentum of the CMG rotor 


H 


w 


all 


b33 


c 

J» 

c 

J' 


c33 


[ J p ] 


K 1 

K 

m 

[K.] 

*P 

Si 

S 

K si 

L 

"P 

m s 

P 

p! 

j 


CMG wheel momentum vector 

Inertia about axis of rotation 

Inertia of inner glmbal and motor rotor about 

gimbal X axis 

CMG base inertia about the Z axis 

Outer gimbal inertia 

Equivalent outer gimbal inertia 

Equivalent outer gimbal inertia about the outer 
gimbal Z axis 

Moment of inertia of the CMG rotor about the 
diameter (X or Z axis) 

Inertia tensor of the pallet 

Moment of inertia of the CMG rotor about the 
spin (Y) axis 

Inertia tensor of the orbiter 
Inertia tensor in vehicle coordinates 
Root locus gain 
Gain margin 

Pallet isolator spring constant matrix 

Control law position gain 

Control law position integral gain 

Control law rate gain 

Shockmount spring stiffness 

Lagrangian adjoint 
Mass of the pallet 

Mass of the orbiter 

Performance index 

CMG output coordinate system (j=l,2,3,4) 


r 

o 

R 

o 

S 

S 

g 

T 


Components of the local vertical unit vector 

resolved in vehicle coordinates 

Vector locating pallet CM to orbiter CM 

Laplace operator 

Friction level at which switching from static 
to running friction occurs 
Torque output from the CMG 
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Time to accelerate to maximum vehicle rate 
for GG maneuver 

Torque exerted on the inner gimbal by the 
outer gimbal 

Maneuver time during constant rate period 
during GG maneuver 

Torque exerted on CMG base by shockmount 
expressed in B-space 
CMG command torque 

Torque exerted on the outer gimbal by the 
CMG base 
Commanded torque 

Total control torque exerted on the pallet by 
the CMG cluster 
Disturbance torque 

Output torque of DCMG 

Disturbance torque acting on the orbiter 
Error torque command 

External torques acting on the inner gimbal 

Friction torque on orbiter 

Friction torque 

CMG gear train torque output 

Gravity gradient torque 

Upper bound on gravity gradient torque 

Gravity Gradient torque acting on the orbiter 
Total maneuver time 
Torque about rotor Y axis 

Torque exerted on the pallet by the shockmount 

Torque for the jth CMG expressed in p* coor- 
dinates 

Running level of bearing friction 


Output torque of SCMG 


+ U tu 


on 


the 


T 

sg 



->* 

T 

w 


T 

wx 

w~space 

X I ,Y I ,Z I 

VV Z P 

W z s 

XIOP 

XPOP 






6 . 
rax 

AI 

A(s) 


e 

g 


-*■ 



Static level of bearing friction 

Sum of the suspension torques acting 
pallet 

Sura of the suspension torques acting on the 
orbiter 

Torque exerted on the w-frame written in w- 
space 

Wheel reaction torque about inner gimbal X 
axis 

CMG rotor (wheel) coordinate frame 
Inertial coordinates 

Pallet coordinates 

Orbiter coordinates 

X vehicle principle axis in the orbit plane 
X vehicle principle axis perpendicular to the 
orbit plane 

Locates pallet end of the ith isolator 

Locates orbiter end of the ith isolator 

Small angle rotation of CMG rotor wrt A-space 
(1=1,3) 

CMG gimbal angle (1=1, inner; i=3, outer) 

CMG gimbal rate (i=l, inner; i=3, outer) 

CMG gimbal rate vector (1=1, inner; i=3, 
outer) (j=l,2, 3,4 for four CMGs) 

Inner gimbal rate command 

Outer gimbal rate command 

Rate of CMG motor rotor (i=l,3) 

Difference in inertia terms 
Control loop characteristic equation 

Pointing offset 

Mass offset of the pallet wrt to the pallet 
initial position expressed in I coordinates 
Threshold level at which friction model 
switches from running to stiction state 

Small angle rotation about the vehicle axes 
Damping factor (i=l,2,n) 


12-5 


Sensed vehicle position 
Desired vehicle position 

Pallet position 

Orbiter position 

Small angle rotation of CMG base wrt the 
pallet 

Lagrange multiplier 

Location of pallet CM wrt orbiter CM 
Initial card on p p , i.e., value of p_ when 
suspension is relaxed 

Relates pallet coordinate system to origin 
of I coordinate system 

Initial condition on p , i.e., value of p 

b S 

when suspension is relaxed 

Relates orbiter coordinate system to the 
origin of the I coordinate system 
Phase margin 

Maneuver angle for desaturation maneuvers 
Sensed vehicle rate 

Acceleration of the ith coordinate frame 
Inertial rate of the CMG inner gimbal 
Inertial rate of the CMG base 

System bandwidth 

Inertial rate of the CMG outer gimbal 
Desired vehicle rate 
CMG gimbal rate 

Compensator undamped natural frequency (i=l, 
Orbital rate 

Spin rate of the CMG rotor 

Inertial rate of the CMG w-frame 
Vector dot (inner) product 


X 

sec 


Vector cross product 
Arc second of rotation 


12.4 T ransformations 

t a T c J Transformation from C-space to A-space 

^Aj T Cj^ Transformation from outer gimbal to inner 

gimbal coordinates 

fb T p'^ Transformation from p’-space to B-space 


w 

Transformation 


Transformation 

coordinates 

w 

Transformation 

dinates 

c x T s^ 

Transformat ion 
coordinates 

tp T v ] 

Transformation 

coordinates 

^•p T p* ] 

Transformation 

p j 

sys tem 

I T ] 
w a J 

Transformation 


from B-space to C-space 
from p' to CMG outer gimbal 

from pallet to inertial coor- 

from orbiter to inertial 

from vehicle to principle axis 

from the p* to the p coordinate 

from A-space to w-space 
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